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(D . Abstract 

£' 

In this paper we study the integrabihty of a family of models with U{1) x SU{N) sym- 
metry. They admit fermionic and bosonic formulations related through bosonization and 
subsequent T-duality. The fermionic theory is just the CP^~^ sigma model coupled to a 
self-interacting massless fermion, while the bosonic one defines a one-parameter deforma- 
tion of the 0{2N) sigma model. For N = 2 the latter model is equivalent to the integrable 
deformation of the 0(4) sigma model discovered by Wiegmann. At higher values of N 
we find that integrability is more sporadic and requires a fine-tuning of the parameters of 
the theory. A special case of our study is the N = 4 model, which was found to describe 
the AdS^ X CP^ string theory in the Alday-Maldacena decoupling limit. In this case we 
propose a set of asymptotic Bethe ansatz equations for the energy spectrum. 
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1 Introduction 



Two-dimensional theories are physically interesting. They are easier to work with than 
the four-dimensional ones yet they exhibit phenomenons like confinement, charge screen- 
ing, dynamical transmutation, etc. A prominent example is the so-called CP^~^ sigma 
model [H El E] that mimics and simplifies many of the expected features of four-dimensional 
confining gauge theories. Microscopically the model is defined by the Lagrangian 

C = K{d^-tA^)z{d'' + tA^)z, (1) 

where z = {zi, . . . , Zn) is an SU{N) multiplet of complex bosons subject to the constraint 
zz = 1 and is an abelian gauge field. It has U{1) gauge symmetry — )■ — d^u, z — )■ 
e^'^z, but no kinetic term for the gauge field. The latter may be integrated out and replaced 
by a local interaction = izd^j^z. The model ([1]) has a single dimensionless coupling 
~ 1/k, it is known to be renormalizable ^ and is asympotically free. To the leading 
order ~ /3olog(/i/A) with A being the dynamical scale and /3o = N/2t[. Similarly to 
the 0{2N) model [SIE], the model ([T]) is solvable at large iV, see [21I3|, and the fundamental 
excitations acquire a mass ~ A. What makes it different, however, is that the fundamental 
excitations are confined by the long-range Coulomb interaction induced non-perturbatively 
by the gauge field. The spectrum has a gap and is populated by "mesons" falling into 
representations of SU{N)/Zn, with the center of SU{N). 

Much less is known about the CP^^^ model beyond the large N limit. Classically 
the model has infinitely many conserved currents and is completely integrable pQ. Un- 
fortunately, these conservation laws are spoiled by anomalies at the quantum level [8] 
and the integrability is not preserved by the quantization. A noticeable exception occurs 
for N = 2 when the theory is equivalent to the 0(3) sigma model [2], [3] and hence inte- 
grable P, Uni [m 112] ■ a more promising class of models from the perspective of integrability 
is obtained by minimally coupling a massless Dirac fermion to the gauge theory ([1]). The 
corresponding Lagrangian is 

£ = K{d, - tA^)z{d^ + tA^)z + t^rid, - ikA^)^ - ^ (^7/.^)' , (2) 

where the fermion has charge k and a self-interaction controlled by the Thirring coupling 
A. As opposed to ([T]) the theory ([2]) has an additional U{1) symmetry associated to 
the conservation of the number of fermions. At the classical level it also has an axial 
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f/(l) symmetry coming from chiral rotations of the fermion. This symmetry becomes 
anomalous at the quantum level and is broken down to a discrete Z2fe subgroup, which in 
turn is spontaneously broken down to Z2 [3]. 

From the viewpoint of perturbation theory both models ([T]) and look similar. For 
instance, to leading order at weak coupling the running of the coupling I/k is oblivious to 
the presence of the fermion and the theory ([2j) is equipped with the dynamical scale A of 
theory ([T]). The most striking difference between the models ([T]) and is non-perturbative. 
As evidenced by the large N analysis [3l |T3] the U{1) gauge symmetry is spontaneously 
broken in (|2]). The phenomenon is similar to the one observed in the Schwinger model; 
the massless fermion is eaten up by the gauge field which in turn acquires a mass of order 
~ fcA/ iyiV(r^rX77r) for large N and becomes dynamical. This leads to the screening 
of the Coulomb interaction at long distance and the fundamental excitations, henceforth 
dubbed spinous, are liberated. A natural question in these new circumstances is whether 
the system can now become integrable. 

In this paper we will give a positive answer to this question. We will argue that for 
any value of and the fermion charge k there exists one value of the Thirring coupling, 
or more accurately one renormalization group (RG) trajectory, for which the system ([2]) is 
integrable and described by the minimal reflectionless U{N) S-matrix. The latter S-matrix 
corresponds to the class II solution in the classification of U{N) invariant factorized S- 
matrix established by Berg et al. 

We should stress at this point that our discussion is not unrelated to earlier consider- 
ations. It has been known for a long time that the integrability of the quantum CP^~^ 
model can be restored by adding massless Dirac fermions [Ul [IS]. Such a model was pro- 
posed by Koberle and Kurak [12] who introduced A^ Dirac fermions with charge k = 1. 
The coupling to the bosonic degrees of freedom was chosen to be minimal and the Thirring 
coupling was set to zero. To verify their proposal the authors of [16] constructed the S- 
matrix to leading order at large A^ and observed it to agree with the large A^ expansion 
of the minimal reflectionless S-matrix of [H]. Further support for the integrability of this 
model was given in [T71 [18] . 

The model which we study in this article is not essentially different from the one 
considered by Koberle and Kurak in [16]. It has however couple of advantages over the 
latter. First of all, theory (|2]) is minimal because it contains only one massless Dirac 
fermion. It is known that only one fermion is effectively active even if multiple massless 
fermions are minimally coupled to a f/(l) gauge field in two dimensions [IH]. This was also 
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observed by the authors of [TB] . From this point of view the model appears as the 
most generic effective theory. Furthermore, the theory ([2]) is renormahzable, at one-loop 
at least, as opposed to models without Thirring coupling. We shall prove that the Thirring 
coupling A in ([2]) is not exactly marginal and starts running at the one-loop level 

A = Aoo-^ + 0(l/«:^). (3) 

The parameter Aqo denotes the UV value of the Thirring coupling A. It is physical and 
may be chosen arbitrarily. Consistency with the factorized S-matrix requires the Thirring 
coupling to run and the parameter Aoo to be fine-tuned. We shall show this by studying 
the free energy density of a gas of spinous at large chemical potential. The computation 
may be done either by using the conjectured exact S-matrix or with help of the standard 
perturbation theory. Matching both calculations is very sensitive to the renormalization 
properties of the Thirring coupling and confirms (|3]). The theory (|2]) may thus be seen to 
provide a refined version of the model introduced by Koberle and Kurak such that it is 
consistent with the usual requirements of a perturbative QFT. 

Another motivation for studying the theory ([2]) is tied to the AdS / CFT correspondence. 
A model of this class has been recently found by Bykov [20j to describe the truncation of the 
AdS4 xCP^ type IIA super-string sigma in the Alday-Maldacena decoupling limit [21] . The 
latter is understood as the low-energy limit of the string sigma model in the background of 
a long string rotating in AdS-^ C AdS^ [23. The effective Lagrangian reads explicitly [20] 

C = K{d^ - tA^)z{d^ + %A^)z + #7^(9^ - 2%A;)i) + ^ (V^7mV^)' , (4) 

with the bosonic variables describing CP^ and the massless Dirac fermion being the only 
extra remnant of the superstring coordinates in the massless limit [23j . The effective model 
dH) belongs to the class ([2]) with = 4, = 2, and Aoo = 0. 

It is of particular interest to study the integrability of this effective model. The sigma 
model on AdS^ x CP^, analogously to the AdS^ x S'^ case, is believed to be quantum 
integrable. We refer the reader to [23] for a recent review. The amount of evidence 
supporting this supposition is significantly smaller than for the AdS^ x S"^ theory. In this 
paper we will give a non-perturbative argument in favor of the integrability of the jull string 
sigma model by showing that the model (jl]) belongs to the integrable subclass of ([2|). We 
shall further propose a set of asymptotic Bethe ansatz (ABA) equations for the energy levels 
of the theory ^ in finite volume, i.e., defined on a cylinder of length L, and discuss subtle 
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effects related to the restoration of the ^ Z2fe/Z2 chiral symmetry. We will see that 
these finite-volume effects require a twisting of the ABA equations. Eventually, it should 
be possible to compare our results with the predictions coming from the conjectured all- 
loop ABA equations for the Aharony-Bergman-Jafferis-Maldacena (AB JM) theory 
[27] . This comparison is beyond the scope of our analysis but would provide a stringent 
all-order test of this conjecture. 

Last but not least, we shall show that the fermionic model ([2]) is equivalent to 

£ = ^(5, - ^AM^' + ^Anz + l^y , (5) 

which defines a one-parameter family of marginally relevant deformations of the 0{2N) 
sigma model preserving the f/(l) x SU (N) symmetry. First indications of this duality were 
given in |2S] where the fermionic model considered by Koberle and Kurak was related to 
the bosonic model ([5]) at the level of their large N effective actions. In this article we shall 
establish the equivalence between ([2]) and ([5j) by means of bosonization and subsequent 
T-dualization. This will lead to the following identification of the parameters 

while the charge k will be found to parametrize a quotient of the model Geometri- 
cally the deformation of the sphere S*^^"^ described by the sigma model ([5]) is also known 
as a squashing, with the round sphere being recovered for = R^. Notice also that in 
the bosonic formulation there is no gauge symmetry, due to the "mass term" ~ A^j^A^^. 
For this reason the model is sometimes termed the "massive" CP^^^ model [29j. It is 
nevertheless more suitable to view this model as a deformation of the 0{2N) sigma model 
since it entails 2A^ — 1 massless fields in the UV, that is one more than there are available 
in the case of the CP^"^ model. The latter is recovered in the limit ~ when one 
degree of freedom decouples. 

In the special case of = 2, the sigma model is equivalent to the one parameter 
deformation of the 0(4) sigma model introduced by Wiegmann and Polyakov [201 EI]- This 
model is well-studied and was found to be integrable both at the classical and quantum 
level [HOlEaiMllMlESllMlEZllMlEg]. Since it belongs to the family of models considered 
here, it offers valuable insights into the non-perturbative physics at smallest non-trivial 
value of A^. This theory remains nonetheless exceptional from the view point of integra- 
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bility. Most of the arguments regarding the generic case put forward in this paper hold 
for N > 2 and for each of them we will point out what makes the case = 2 so special. 

This paper is structured as follows. In Section |2] we present the model in its two 
formulations and relate them via bosonization and T-duality. We also discuss their renor- 
malizability and infrared physics. The classical and quantum integrability of the model 
are subject of Section |3l In Section S] we perform a direct test of the proposed S-matrix 
at finite values of N by computing the free energy density of the model at finite chemical 
potential. In Section |5] we apply these results to the particular case relevant to string the- 
ory and construct the associated ABA equations. Section |6] contains concluding remarks. 
Technical details of our analysis are deferred to the various appendices. 



2 The models 

In this section we introduce the fermionic and bosonic formulations, discuss their renormal- 
izability and establish their equivalence. We start with the bosonic model whose physics 
may be more easily exposed. 



2.1 The bosonic model 

The bosonic model is defined by the Lagrangian ([5]), with the A^ complex scalar fields 
z = {zi,...,zn) spanning the unit sphere S^'^~^, i.e., satisfying zz = 1. An equivalent 
formulation is obtained if one decides to integrate out the dummy field A^. One finds 

C = ^id^- iB^)z{d^ + iB^)z + "j-B^B^ , (7) 
47r An 

where B^ = izd^z. The first term in ([7]) has a local U{1) symmetry and coincides with the 
Lagrangian of the CP'^'^ model. The last term breaks this gauge symmetry down to the 
global U{1) subgroup and geometrically describes a circle of radius r fibered over CP'^~^. 
For generic values of r the model ([7]) has thus U{1) x SU{N) symmetry. The symmetry 
is enhanced if the two radii are equal r = R when the model becomes equivalent to the 
0{2N) sigma model. In the limit r — )■ the small circle shrinks to a point and the theory 
reduces to the CP^~^ sigma model. The regime considered in this paper is bounded by 
these two special values ^ r ^ R. 

The Lagrangian ([7]) makes the geometry of the problem manifest and will be the start- 
ing point for performing the T-duality in the next subsection. It is nonetheless entirely 
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equivalent to 

£ = nd^zd'^z + Kr]{zd^zf , (8) 

where r] and k are dimensionless parameters, related to the radii by = Ann , = 
Attk{1 — 1]). This form makes the relation to the 0{2N) model evident and is sometimes 
more convenient for semiclassical computation. Notice that for the afore-stated domain of 
R the parameter rj takes values in the interval ^ rj ^ 1. The boundary values rj = 0,1 
correspond to the 0{2N) and the CP^~^ model, respectively. 

As in the 0{2N) sigma model, the radii start running at the quantum level. To reveal 
this effect one should proceed with the renormalization of the model ([7]). This analysis can 
be done perturbatively as an expansion at large R and was performed in detail in [29l [33] . 
The theory was found to be asymptotically free in the domain ^ r ^ R and the RG 
equations were constructed up to two loops. They take the following form 

for R^ 1 and fixed r. Both radii are monotonically increasing functions of the renormal- 
ization scale fi, but while R grows without any bound, the radius r reaches a finite value in 
the UV. Put differently, in the (r, R) plane these equations describe a flow of trajectories 
ending on a line of fixed points at = oo, see Figure [U This line is parameterized by 
the UV value = r(oo) of the radius r, which by definition is RG-invariant. The second 
RG-invariant is the dynamical scale A which controls the running of the radii. Explicitly, 
the solution to ([9]) is 

R^ = 2N\og (/i/A) + 2(2 - p) log log (/x/A) + o(l) , (10) 

and 

^- ^ ^^'^ {N-l){2-p) 

- 2A^ + W NR- + ^^^/^ ) • ^^^^ 

Here we have introduced the parameter 

P=— , (12) 
^ 2iV' ^ ' 

which will turn out to be convenient when discussing the physical properties of the model. 
Please note that the point p = corresponds to the CP^~^ sigma model. The 0{2N) 
sigma model, on the other hand, is formally recovered when p — )■ oo. 
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1/R 



Figure 1: A specimen of the solutions to the RG equations in the asymptotically free 
domain ^ r ^ R. The upper line correponds to the 0{2N) trajectory r = R and 
the lower one to CP'^^^ with r = 0. The arrows point toward increasing values of the 
renormalization scale fi. All trajectories in the band ^ r < R end on the line 1/R = 
at = oo. 

What can one infer about the non-perturbative infrared features of the model? A 
standard method to derive hints about the physics beyond the perturbation theory is the 
large N expansion. This study has been conducted for the model at hand in and 
closely parallels the one for the dual fermionic model [3l [13]. For the sake of the large N 
analysis, it is convenient to introduce two auxiliary fields: the gauge field as in ([5]) and 
a scalar field a which implements the constraint zz = 1. As in the case of the 0{2N) sigma 
model, it is found that a acquires a vacuum expectation value giving a mass m = A to the 
2N particles. These are our spinous and antispinons transforming in the fundamental and 
antifundamental representation of SU{N), respectively. They have opposite charge with 
respect to the U{1) group. The interactions, which are weak at large A^, are mediated by 
the gauge and the scalar field. The vector interactions play a prominent role in the physics 
of the model and distinguish this theory from the 0{2N) sigma model. They are controlled 
by the gauge field propagator 



N{A{k^) -p) 



—ITT 




+ 



Npk^ ' 



(13) 
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with the function ^(/c^) defined for A;^ < by 



1 I 47772 / iU2 \ U2 

Me) = ^arccc«h(l -_)=_ + OikVrn') . (14) 

and otherwise by analytic continuation. Its structure is reminiscent of the Proca propagator 
for a massive gauge field with mass term controlled by the parameter p. For p = 0, when the 
model is equivalent to CP^~^, it has a pole at A;^ = leading to a linear confining potential. 
The spectrum consists of mesons, which are spinon-antispinon bound states |2l |3]. For 
< p < 1 there is a pole in the propagator (I13p at finite value of k"^. The force is short- 
ranged and the spinous are liberated. As long as p is very small in value, the mass of 
these particles is still much larger than that of the gauge field, which can be viewed as a 
singlet bound state. A transition occurs at p = 1 when the lightest bound states, including 
the gauge field, hit the two-particle threshold. This point is special from the integrability 
perspective, as already pointed out in [28] building on [16j, and we shall come back to 
it later on. Finally, the range p ^ I corresponds to repulsive interactions and no bound 
states. We refer the reader to [28j for a more detailed discussion of the large N theory. 

2.2 T-duality 

To unveil the relation between the bosonic and the fermionic theory we need to T-dualize 
the former along its fibered circle^ To perform this transformation it is easier to work with 
the Lagrangian ([7]) and use local coordinates which trivialize the U{1) bundle. One can 
write for instance Zi = e^"^ z[ where ■(9 is a compact boson with period 27r, {} ~ -(9-1-27?, and z[ 
is a set of local coordinates on CP^~^ . In these coordinates the one- form i?^ decomposes 
into vertical and horizontal parts, 

B^ = -d^^^h^, (15) 

where is a local one-form on CP^^^. In this formulation the U{1) symmetry is imple- 
mented as translations of d and this isometry allows us to T-dualize along this direction. 
Notice in this regard that all the dependence on d resides in the second term in ([7]) and is 
manifest in f llSp . To perform the T-duality we will therefore focus only on this part of the 
Lagrangian, which we denote as C fiber- 

To accommodate for the fermion charge fc, one actually needs to consider a quotient 

^We are most grateful to D. Andriot for helpul discussions on this topic. 
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of the bosonic model obtained by identifying 



z ~ e^'^l^z . (16) 

It is equivalent to redefining the field d d/k while keeping its 27r-periodicity. Taking 
this into account we need to dualize 

2 

^Me. = £^(9^^-A;6^)^ (17) 

with h^j, treated as a background field. A detailed discussion of its T-duahzation can be 
found in [lO]. Further references on the subject with direct connection to the model under 
consideration are [HI The outcome is simply 

Cf^ber = %d^^d^^ + ^ e^^d.B, , (18) 

now written in terms of the T-dual compact field ip. In the above expression we have made 
use of the following relation e^^d^h^ = e^'^d^By with e^^ being the antisymmetric unit 
tensor. Note also that in our notations the field (p has period a/vt, which is a convenient 
choice for the mapping with the bosonized fermion. The charge k could be eliminated by 
redefining the field ip ^ ip/k, hence changing the period of p. It is associated to a discrete 
Zfc symmetry which is spontaneously broken [2j and is irrelevant in infinite volume. It leads 
however to interesting effects when the model is put on a finite cylinder, as we shall see in 
Section O In view of this fact and in order to facilitate the mapping with the fermionic 
model we prefer to keep the charge k explicit as in 0181) . 

As expected, the coordinate transformation leading from (|T7|) to f|T8l) is non-local. 
Explicitly it reads 



<^(a,r) = ^ j da'Ma',r), (19) 

where = r'^B^/ {2ti) = r^{kb^ — d^i!}) / {2nk) is the U{1) current. An equivalent local 
form is given by 

J^ = -^e^,9>. (20) 



The U{1) current is then automatically conserved in this picture and the associated charge 
is topological for the p> field. A further interesting observation is that the T-dual bun- 
dle is trivial, as opposed to the original one. The topological structure of the orig- 
inal bundle translates however into interaction between p> and the topological density 
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e'"'d^,B^ = it^^'d^zd^z of the CP^-^ model. 

Eventually, we can reinstate the dummy gauge field and write the T-dual of the 
model ([5]) as 

C = —D^-zD^z + ——^k^d.^d^^ + -^e^'d.A, , (21) 

where is the covariant derivative built with help of A^. This is the form we will need 
to make contact with the bosonized version of the fermionic model. It is amusing to see 
how the 0{2N) model is recovered in this picture. The latter corresponds to r = R and 
at this point the field (f in ( 12T|) has no kinetic term. It becomes a Lagrange multiplier 
implementing 

e^^a^A, = 0. (22) 

This in turn means that is pure gauge and can be eliminated by a local rotation of 
the z^s. Hence, locally, we recover the Lagrangian of the 0{2N) model. A refined analysis 
would reveal that it is actually a model with target space S^^^^/Zk. In the opposite limit 
r — )• the boson if decouples and we are left with the CP^^^ model. 



2.3 The fermionic model 

The fermionic extension of the CP^^^ model is defined by the Lagrangian ([2]). Apart from 
the inherited coupling k, it contains two additional parameters: the fermion charge k and 
the Thirring coupling A. They are both dimensionless. In this section we will relate them 
to the parameters of the bosonic theory, i.e., we will derive the equations ([6]) quoted in the 
introduction. 

To make contact with the previous formulation we need to bosonize the fermion [12] . 
The bosonized form of the model without Thirring coupling was already given in [5]. Using 
the standard identity [12] 

= ^e^y&'ip , (23) 

where stands for the bosonized fermion, one sees immediately that the net effect of the 
Thirring coupling is to redefine the kinetic term of the bosonized theory. Explicitly, the 
disguised fermionic Lagrangian takes the form 

Hp = ^^^d,^d^^ + ^e'^^d.A, . (24) 

2 y/TT 

It is identical to the T-dual version of the quotient of the bosonic model fl2T]) if the 



12 




Figure 2: The relation between the models. 

couplings are identified as in ([6]). The relation between the bosonic and fermionic models 
is summarized in Figure |2i 

The relation established above allows us to translate the physics of the bosonic model 
into the fermionic language. Firstly, the Thirring coupling A, being essentially the same as 
the radius r of the bosonic model, is not exactly marginal and runs with the renormalization 
scale. In the UV it approaches a constant value Aqo which is a RG-invariant parameter of 
the model. Secondly, we learn that the true deformation parameter of the fermionic model 
is neither k nor X^o on their own, but the combination 

^ = iV(l + Aoo/vr) • ^^^^ 

This quantity has to be identified with the parameter p = r'^/2N introduced for the 
bosonic model. The relevance of this parameter for the fermionic model could be inferred 
from the large N analysis by generalizing the argumentation of [3] to non-zero value of the 
Thirring coupling. It is important to stress that the renormalization of this interaction is 
subleading at large and may be ignored to leading order. One finds that p controls the 
mass of the gauge field, which is the only physical parameter in the theory. Varying the 
fermion charge k (assumed to be quantized) at a given value oip is equivalent to performing 
a Zfc-quotient of the bosonic model. As mentioned before, this operation has no effect on 
the model in infinite volume and, for instance, does not affect integrability. 

This is the picture we inherit from the duality with the bosonic model. In what follows 
we will argue that one also recovers it directly from the fermionic formulation. 
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The fermionic model has naively more symmetries than the bosonic one. It has U{1) x 
SU{N) symmetry from rotations of the z's and U{1) x U{1) from vectorxaxial rotations 
of the fermion. Due to the U{1) gauge symmetry, one U{1) subgroup is removed. This 
may be observed directly at the level of the Gauss law 

2iKzD^z = k-ijj'y^'ilj , (26) 

where = 2inzD^z is the bosonic U{1) current. There remains only one extra sjTiimetry, 
which has no equivalent in the bosonic formulation. It derives from the axial U{1) rotation 
of the fermion, 

^ ^ e^"^^?/' , (27) 

with 75 = 7o7i. Fortunately, it is purely classical and is spoiled by the anomaly upon 
quantization. This is manifest in the bosonized form of the model, where the axial U{1) 
transformation fl27|) translates to ip if + a/ a/tt and is not a symmetry of fl2Tl) . What 
remains is a discrete chiral subgroup Z2fc which is spontaneously broken down to Z2 |3]. 
They are generated hj a = n/k and a = ir, respectively. Up to the Z2 symmetry, the 
quantum symmetries of the fermionic model are then identical to those of the bosonic 
model quotiented by = Z2fc/Z2. We recall that discrete symmetries play no role in 
infinite volume. So before we put the theories on a finite cylinder in Section O the reader 
might simply ignore them. 

The second point we wish to understand directly from the fermionic picture concerns 
the renormalization of the Thirring coupling A. In the absence of interactions between 
bosons and fermion, i.e., when k = 0, the Thirring coupling is exactly marginal. We shall 
prove below that this is not the case anymore if the interactions are switched back on. 

Before we do that we want to comment on the renormalization of the electric charge 
k, or more precisely on the lack thereof. To make this point manifest in the fermionic 
picture one needs to take note of equation ( l26i) . This equation is an identity between the 
bosonic and fermionic U{1) currents. They are physical observables and hence UV finite. 
As a consequence, the charge k, which is just a proportionality factor between these two 
currents, should not renormalize. 

Let us now turn to the Thirring coupling. To reveal its dependence on the renormal- 
ization scale we can consider the two-point function (^ip'~^^ip[x)%lj'~^yip{Q)') . The interesting 
property of this correlator is that it is sensitive to the Thirring coupling at lowest order in 
perturbation theory. It also follows from our previous discussion that it should be finite. 



14 



This requirement is easily seen to be in conflict with a direct perturbative computation 
unless the Thirring coupling gets renormalized. This proves the assertion that A in (|2]) is a 
running coupling. The one-loop analysis is performed in Appendix |X1 From it we deduce 
that the complete RG equations read 

''a^=2^ + 0(l/'^'- = _ + 0(1/.') . (28) 

They agree with the bosonic ones ^ upon the identification The equation for the 
bosonic coupling k is the same as in the CP^~^ model. This will only be true at this 
order and is actually a consequence of the fact that the correlator of fermionic currents is 
finite. Notice also that although the Thirring coupling does not appear explicitly in theses 
one-loop /3-functions, no assumption was made about its actual value in the derivation. 
In other words, the equations (128|) are valid to leading order in the 1/k expansion, but 
for arbitrary value of A. We should nonetheless require A > — vr for consistency [12] with 
A — > — TT corresponding to the 0{2N) limit. Finally, we note that the solution to (^E^ is 
given by 

TV k'^ 
= — log (/./A) + 0(loglog (/i/A)) , A = Aoo - ^ + 0(1/ k') , (29) 

where the parameter Aqo is the UV value of the Thirring coupling. 

At the end it is quite natural to find that the Thirring coupling runs in this theory. The 
associated vertex is fully consistent with the symmetries of the theory. In fact, it is the 
only marginal interaction that we can add to a massless Dirac fermion minimally coupled 
to the CP^"^ model, up to an irrelevant 6 term (see discussion in [^). One could also 
attempt to add mixed or bosonic current-current interactions. Due to the Gauss law fl26l) . 
they would not essentially differ from the Thirring interaction. The argument is somewhat 
heuristic but one can check explicitely that the two possible extra vertices are effectively 
redundant. The model ([2]) is hence the most general one with the given symmetries and 
thus the most likely to define a renormalizable theory. 



2.4 The case N = 2 

Before we discuss the integrability of the U{1) x SU{N) model it is worthwhile to briefly 
summarize what happens at = 2. This particular case provides a remarkable illustration 
of the previous generic considerations and has been extensively studied in the literature 
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[SniElESlIMlESlESlEZlESlES]. 

To make contact with the standard formulation one just needs to repackage the two 
SU (2) doublets into a single SU{2) group element 



= • (30) 

Two natural sets of transformations on Vt are left and right SU{2) rotations, Q — )■ LQR, 
with L,R,E SU{2). The left SU{2) group rotates the multiplets Zi and eijZj, while the 
right SU{2) acts on {zi, —^2) and {z2, Zi). The right transformations are not all symmetries 
of the model for generic value of rj. This becomes manifest after rewriting the Lagrangian 
in terms of the right SU(2) currents = Q~^dnQ, 

Cn=2 = ~tr j,f + ^ (tr j^asf . (31) 



In the above a, denote Pauli matrices. One easily recognizes in fl3T|) the deformation of the 
SU{2) Principal Chiral Field (PCF) first considered by Wiegmann and Polyakov [301 EI]- 
For 7] = the model is identical to the SU (2) PCF model, well-known to be equivalent to 
the 0(4) sigma model. Away from this point the second term in the Lagrangian breaks 
the right SU{2) symmetry down to U{1). 

The model (131!) is classically [321 [SU [361 EH EH EH] and quantum-mechanically [301 ESI 
EH ES] integrable for any value of the deformation parameter p. Its S-matrix was derived 
by Wiegmann [50] and found to be given by the tensor product 

§ = §5C/(2)®4g, (32) 
where the first factor is the minimal SU{2) S-matrix and the second one coincides with 



the sine-Gordon S-matrix with the parameter The value p = 1 corresponds to the 
free Dirac point for the sine-Gordon S-matrix. At this point the S-matrix ( 132|) simplifies 
considerably and becomes the minimal reflectionless U{2) S-matrix. For p < 1 the theory 
has spinon-antispinon bound states [501 E3]- These breathers come in pair (singlet, vector) 
w.r.t. SU{2) and are degenerate in mass. The lightest breathers consist of a C-odd singlet, 
which can be seen as the dual field if or equivalently e^i^d^zd'^ z, and a C-even vector created 
by Xi ~ Z(7iZ. At small enough p these bound states are the lightest excitations in the 



^We refer the reader to [43J for a definition of this parameter in the context of the sine-Gordon theory. 
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spectrum. At p = they decouple from one another, with the singlet becoming a free 
massive field and the vector describing the 0(3) CP^ sigma model. Up to the accidental 
degeneracy, this picture is remarkably close to what is expected to happen for higher 
values of A^. In the latter case, however, integrability will place severe restrictions on the 
deformation parameter p, as will argue in the following section. 

Surprisingly, the S-matrix (132|) cannot be found in the U{N) classification of Berg et 
al. [11]. A natural question is whether this classification could be missing a one-parameter 
family of U{N) solutions. Such a family would be particularly significant for the class of 
models we are considering. We clarify this point in Appendix [Bl where it is shown that the 
Wiegmann's solution does not admit a lift beyond A^ = 2 in agreement with the conclusion 
of Berg et al. [H]. We nevertheless found a caveat in their classification for N = 2, which 
allows us to embed the S-matrix ( 132|) into the space of A^ = 2 solutions. 

It is also worthwhile considering the fermionic formulation of the N = 2 model at the 
special point p = 2. To understand what is special about this value, we recall [37] that 
the complete symmetry of the model is SU{2) x SUq{2), where SUq{2) is the quantum 
deformation of U{su2) with q = exp(z7r/p). This symmetry is somewhat visible in the 
decomposition ( 15^ if one recalls that SUq{2) is the quantum symmetry group of the sine- 
Gordon theory [H]. It is however not manifest at the level of the Lagrangian and for q ^ 1 
it is implemented by non-local currents [37], see also 03]. At p = 2 the SUq{2) symmetry 
becomes isomorphic to the centrally extended M = 2 SUSY algebra with 4 supercharges. 
Interestingly enough, this SUSY algebra can be realized linearly in the fermionic model. 
This requires k = 2, when the theory has two conserved spin 3/2 supercurrents proportional 
to eyZ^D^z^'^'^'y^ip and its complex conjugate. What is this model? It turns out that this 
is the SUSY CP^ model. This is not the standard formulation of this theory, which 
by definition is equipped with a SU (2) doublet of Dirac fermions constrained by ZiXi = 
ZjXj = 0, nevertheless it is equivalent [15]. The key observation is that one can solve these 
constraints in terms of a single Dirac field ip defined such that 

Xi = i^eijZj , Xi = i'f^ijZj . (33) 

This fermion has charge = 2 in our notations. The supercurrents mentioned above are 
then nothing else than the supercurrents of the SUSY CP^ model expressed in terms of the 
fermion ip. Interestingly, the Thirring coupling for this model is fixed by supersymmetry 

A = --. (34) 
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As a consistency check of our previous results, we observe that it is in perfect agreement 
with ( 125|1 and f l28|) ii N = p = k = 2. This concludes our discussion of the case N = 2. 

3 Integrability 

In this section we study integrability of the fermionic and the bosonic model for generic 
N. We will find that it is a sporadic phenomenon requiring fixing the parameter p. 

3.1 Classical integrability 

We will start with the classical integrability. It turns out that only the fermionic model 
is classically integrable for generic value of the deformation parameter. This may appear 
contradictory at first, as we have just shown the equivalence of the both models. We will 
resolve this paradox in Subsection 13.1.21 

3.1.1 Bosonic model 

Since the actual analysis is slightly technical, we will first offer a summary of the results. 
Customarily, classical integrability of a physical system relies on the existence of a non- 
abelian current j^^ , which is both flat and conserved, 

d^i^ - d^^ + []„ I] '' = 0, d^j^ = . (35) 

Its existence allows one to construct the Lax connection 

^/^(^^ = ^ ' ^^^^ 

which is automatically flat for any value of the spectral parameter x. The latter property 
guarantees the conservation of infinitely many non-abelian charges, which are manifestation 
of the integrability of the model. They may be generated with help of the monodromy 
matrix defined as the path-ordered exponential of the Lax connection 

M(x) = Pexpj daL^ix) . (37) 

As will be shown below, for the bosonic model the current may only be constructed in 
three exceptional cases. These special cases correspond to fixing one of the two parameters 
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of the classical theory: rj and A^. The case of 77 = is perhaps the most obvious one since 
the Lagrangian ([5]) reduces to that of the 0{2N) sigma model, which is well known to be 
classically integrable. Due to the symmetry enhancement this case is actually not covered 
by our analysis for generic A^. The value rj = 1 corresponds to the CP^^^ model and is 
thus also an integrable case. Our analysis below shows that as long as is arbitrary, there 
are no other special values of rj. The situation is very different if and only if = 2. In 
this case the classical integrability is present irrespectively of the value of the parameter rj. 
This is due to a special kinematical relation which ceases to hold for higher values of A^. 

To substantiate the above picture, we will construct the most general conserved current 
and impose the flatness condition. For an easier comparison with the fermionic theory, it 
is convenient to work with the formulation of the bosonic model and to introduce the 
couplings 7] and k defined previously. The equations of motion are then identical to the 
ones describing the CP^^^ model 

D^D^'z' + D^zD^'z 2' = . (38) 

The sole difference is that the gauge field 

A^, = irjzdf^z (39) 

now depends on the deformation parameter rj. Since we are discussing the classical theory, 
we can further set k = 1. The currents corresponding to the global U{N) symmetry of the 
model can then be written as 

j7 = zW^z^ - D.z'z^ . (40) 

Notice that these currents are not real because we rescaled them by a factor —i. It is also 
convenient to single out the U{1) component 

j'^ = ti3^ = 2zD^z, (41) 

which was previously denoted — z in Section [2l An important remark regarding this 
current is that it is not axially conserved 

e^^S^j; = 2(1 - r^),^'^D,-zD,z = 2(1 - r/)e^^5^z9,z . (42) 
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An exceptional case is rj = 1 when the U{1) current is actually zero. 

The most general conserved current is a combination of the above currents and a topo- 
logical term [311 ES] 

Ji:^ = aji^ + bj;^^^+e,,d^r, (43) 

where a, b are both constant and / is a matrix function. For 77 = the symmetry is enhanced 
to 0{2N) and hence the theory has extra conserved currents. We shall not consider this 
case in what follows. It is important to point out that for rj = 1 the parameter b is irrelevant 
since j'^ = 0. Finally, the function / in fj^S]) should have mass dimension zero. Its most 
general form is then given by 

P = f{zz)z'z^ = cz'z\ (44) 

where c = f{zz) = /(I) is arbitrary. If one can choose a,b,c in such a way that the 
current fH3|) is flat, it will immediately allow us to construct the Lax connection fl36|) . 

Let us check whether such a choice is possible. Using (138|) . the current ( 143|) is easily 
found to obey 

4 

d^i^ - + [j„ I] = e,. ^^^k - ^^.(1 - a)nr ■ (45) 



k=l 



Here the operators O^^ are defined by 

= e^^'Dpz'D^z^ , O'J = eP^DpzD^z fz? , 

= ef'zDpZ d^iz'z^) , O'J = e^^D^zD^z 6'^ , 

with O3' being identically zero if r/ = 1. The coefficients in f H5]) are given by 

«! = — a(2 — a) , a2 = a{2r] — a) — , 
as = — , a4 = (r] — 1)6 . 



(46) 



(47) 



If the current is flat, then the last term in the r.h.s of (H5|) should vanish separately, as 
it is the only term which is even under the charge conjugation z z. This leads to two 
possibilities: a = 1 or c = 0. Imposing now that the sum in fH5|) vanishes term by term 
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yields two solutions to the flatness equation (153]) : 

a = l, c^ = l, '?? = 1, or 

(48) 

a = 2, c = 0, 7] = 1 . 

These two solutions produce the same charges after expanding the monodromy matrix ( 137|) 
at large x and are consequently equivalent. The solution with a = 2 is the conventional 
expression for the Lax connection of the CP^~^ sigma model. 

As it is clear from ( HHj) no solution is found for rj ^ 1. This disappointing conclusion can 
be avoided if and only if the operators fH6l) are linearly dependent. In fact, the above-listed 
operators form a basis of dimension two operators that are antisymmetric with respect to 
charge conjugation. Any reductiveness of this basis would imply 

4 
n=l 

for some values of c„, n = 1,2,3,4. The complete set of constraints on the coefficients c„ 
may be derived by contracting with 5*-' , z^z\ z^DaZ^ and DaZ^ D°'z\ respectively. One finds 
the following set of equations 

(AT - 2) C2 = , (Oi + O2 + O3 - O4) C2 = . (50) 

Clearly, if any, the non-trivial solution exists only for N = 2. An explicit parameterization 
in this case allows to verify that flSOj) is satisfied for any C2. The current fl43p is then found 
to be fiat for any value of r/ if a = 1, = r/, and h = —1. The two solutions c = ±77 
are related by charge conjugation and generate the same charges. This fiat and conserved 
current, albeit in a different formulation, has been found in [33] and independently in [S7] . 

The above construction suggests that the classical integrability is not present when 
1] ^ 1 and N > 2. Before discussing the fermionic formulation let us delve into the 
possibility of having more general Lax connection than fl36p . We are thus led to study the 
most general solution to the fiatness condition 

d^L, - d,L^ + [L^, L,] = , (51) 

modulo the kinematic constraint zz = 1 and the equations of motion (!38l) . For r] ^ 1 
and A^ > 2 we found that no non-trivial solutions exist, up to the gauge transformation 
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generated by the local unitary matrix = 6^^ + az'^zK This corroborates the claim that 
there is no classical integrability for generic values of rj and A^. 

3.1.2 Fermionic model 

The situation looks differently if one chooses ([2]) instead of as the starting point. We 
have seen that the CP^~^ model is integrable. The Lagrangian ([2j) is an extension of that 
model with a massless Dirac fermion and it is known that one can construct a flat and 
conserved current in this case |17j . Let us briefly present this construction. The equations 
of motion for the fields are the same as before fl38|) and so is the expression for the U{N) 
currents. The gauge field now also accounts for the fermion 



Note that the bosonic U{1) current, being proportional to ip'j^iip, is axially conserved at 
the classical level 



This makes the analysis very similar to the CP model, apart from the fact that the 
f/(l) current does not have to vanish in the fermionic model. Our previous analysis carries 
over to this case, but the coefficients au appearing in ( H5|) now need to be replaced by 
their values f H7|) evaluated at r/ = 1. The immediate conclusion is that the conserved 
current f H5]) is flat for = a = 1 and for any value of h. The latter parameter controls 
the f/(l) part of the mondromy matrix which is not of interest here. Fixing 6 = and 
choosing any one of the two values c = ±1 one easily generates the non-abelian charges 
constucted in [T7], by expanding the mondromy matrix fl57|) at large x. Notice that the 
above construction is independent of the value of the Thirring coupling, since it only entails 
the axial conservation law fl53|l . 

As already pointed out at the beginning of this discussion, the classical integrability 
of the fermionic model does not contradict the non-existence of the corresponding current 
for the bosonic formulation ([5]). The reason is that the bosonization is not an innocent 
classical transformation: it automatically accounts for the axial anomaly, which in turn 
spoils the classical integrability of the bosonic model. 

The above observation has far-reaching consequences. If the classical bosonic theory 
incorporates the effect of the axial anomaly of the fermionic model, should we not conclude 
that its description is more reliable and that integrability can only occur at ?7 = 1? We 



= izd^z - —ip'jf.ip . 



(52) 



e^'-d.jl = . 



(53) 
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believe this not the case. The reason is that the classical bosonic theory only resums a 
subclass of quantum corrections of the fermionic model. To make this point more precise 
we observe that the classical fermionic and bosonic descriptions probe different parts of 
the parameter space. This becomes clear after recalling the relation between the couplings 
of both theories, 

where dots stand for quantum corrections ~ The fermionic model becomes classical 

for K ^ 1. In this limit the parameter r/ takes the classical value t] = 1 for any finite 
value of p. Moving away from this value requires p ~ k. Since we found that classical 
integrability is not a property of the bosonic model for > 2, it is reasonable to infer that 
the fermionic theory is not integrable at the quantum level for arbitrarily large values of 
p. Otherwise the phenomenon should be visible in the corresponding semiclassical regime 
of the bosonic theory. A remarkable illustration is the case N = 2, for which the theory is 
quantum integrable for any p and classically integrable for any rj. 



3.2 Quantum integrability 

We have found that classical integrability of the fermionic model does not imply classical 
integrability for its bosonic dual. Moreover, it seems unlikely that the fermionic model 
will remain quantum integrable at large values of p. In this section we will argue that 
this conclusion should also apply for more generic values of p, hence leaving little room 
for integrability when N > 2. The large N analysis will reveal that integrability is only 
achieved at p = 1 and p = oo. 

3.2.1 Counting and fine tuning 

An elegant way of providing evidence for the quantum integrability of a certain class of 
models offers the counting argument advocated in [46j. The underlying idea is the following. 
Observe that the classical theory is scale- invariant. Then in light-cone coordinates the 

conservation of the stress-energy tensor takes the form 9+T = 0. This for example 

trivially implies 

5+(T_r = 0, (55) 

for any positive integer n. If there exist some axially conserved currents in the theory they 
may also be taken into account in the product. The scale invariance of ([2]) is broken at 
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the quantum level. But for integrable models the above conservation law is expected to 
become deformed not spoiled. Indeed, as long as the deformation is of the form 

a+((T__r + Fi) = a_F2, (56) 

for some operators Fi and F2, it leads to a conserved quantum charge. It was proposed 
in [l6] to compose a list of all quantum anomalies modulo kinematical constraints and 
equations of motion that can appear on the r.h.s of flS^ . The next step is to list all 
admissible operators of the form and 9_ (...). If the both lists match, every anomaly 

may be expressed as a divergence of some operator, consequently allowing to write fl56|) at 
the quantum level. 

We have classified all possible anomalies and divergences for two of the higher conser- 
vation laws and found that in both cases there is one unmatched anomaly if > 2. We 
present the lists for one of the higher conservation laws in Appendix O In principle, there 
might exist a higher charge other than the ones considered by us for which the counting 
would go through. However, in light of the fact that one of the conservation laws we have 
studied was sufficient to argue for integrability for variety of models [16] we consider this 
as a strong evidence against the integrability of ([2]) for generic values of and p. 

We would like to stress that the counting does not depend on the values of the parameter 
p. Since the lists differ only by one entry and the theory has one free parameter p, it may 
happen that the coefficient of the unmatched anomaly cancels out for some adjusted values 
of p. We suspect that this mechanism will restore the integrability for p = 1 and shall give 
evidence for it in the following. It should also be observed that at p = 00, when the 
symmetry is enhanced to 0(2A^), the counting argument leads to correct conclusions [46] . 
Notice finally that the case N = 2 appears again to be exceptional. At this value one 
needs to impose an additional kinematical relation [16] , which reduces by one the number 
of anomalies in the list. The higher conservation law is then guaranteed to exist at the 
quantum level, independently of the value of p. 

3.2.2 Large A^ analysis 

To shed light on the possible values of the parameter p for which the model becomes 
integrable, we shall look at the spinon-(anti)spinon elastic scattering amplitudes at large 
A^. Whether the theory is integrable or not, the two-to-two scattering processes are always 
constrained by the kinematics in two dimensions. The two incoming momenta ought to be 
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conserved separately and the amplitudes are functions of a single Mandelstam invariant. 
In the case of the U{1) x SU{N) model we have three distinct elastic processes to consider. 
They are associated with the S-matrix elements 



{P3k,P4imPu,P2j) = Sf^{d)5{p,~p;)5{p2-P,) + Sf^{0)5{pi-p,)6{p2-p,) 
(p3fc,P4z|Sbii,p2i) = Fg{9)5{pi-P3)5{p2-P4) + Bl';{9)5{pi-Pi)5{p2-P3) 



(57) 



where the asymptotic states are normalized as {p'\p) = S{p — p'). Here p,p are spinon, 
antispinon momenta and 6 = 2 arcosh(s/4m^) is related to the square of the center-of- 
mass energy s. 

The symmetry of the problem permits to decompose the matrices S, F and B into 6 
scalar amplitudes [H] 

F^^{e)=t,{9)5l5] +t2{e)5,,5^\ (58) 

They are not all independent due to the crossing symmetry between s and t = Am? — s 
channel 

Ml 2(m -9)= h^2{0) , n{i7i -9)= r2{9) . (59) 
At large these amplitudes read 

ui{9) = U{9) = 1 , U2{9) = t2{9) = n{9) = r2{9) = , (60) 

as in any weakly coupled (free) theory. The computation of the leading corrections is a 
direct application of the Feynman rules, which are the same as for the CP^^^ model [2] but 
with the gauge- field propagator replaced by (fT3l) . Due to the crossing symmetry (l59l) we 
can further restrict ourselves to the evaluation of ui^2{9) and ri{9). The relevant Feynman 
diagrams are depicted in Fig. [31 One should pay attention to the fact that conventional 
Feynman algebra computes the amplitude A^fj enterinjf 

{p\uM^ - l\pu.P2,) = l^KE - E')6{p - p')M':j{9) , (61) 

where E and p are total energy and momentum, respectively. Hence when relating the 
Feynman amplitude A4^j to the one introduced before, the Jacobian associated to the 



^We recall that we are working with the normalization {p'\p) = S{p ^ p'). 
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Pl,i Pl,k Pl.i P2,l Pl.i P2,k 



V2,j P2,l P2,j Pl,k P2,j Pl.l 

Figure 3: Tree level Feynman diagrams for the scattering amplitudes ui{d), U2{0) and ri(6'), 
respectively. Zigzag lines stand for both scalar and vector exchange. 



change of measure has to be taken into account. It leads to 

Sm = ^'A + ^^^ (62) 



and similar expressions for the spinon-antispinon backward and forward S-matrix: Bfj{6) 
and Flj{9). Summing up the diagrams in Figure [3] gives 

^ in f 2m^ s-2m^\ ^ / 1 

-iG- to + 2(p-l)V^ )+'^(i^ 

In each of these equations, the first /second term in brackets originates from the scalar/vector 
exchange. In the limit p — )■ oo the vector exchange is suppressed and we are left with the 
0{2N) sigma model results written in SU{N) variables. In the opposite limit, when p — )■ 0, 
the spinon-spinon scattering amplitude ui{6) diverges. This phenomenon is a manifestation 
of the charge confinement in the CP^~^ sigma model. 

For generic value of p the scattering amplitudes fl63|l are not consistent with factorized 
scattering, which is the salient property of integrable theories. It can be shown [I^ indeed 
that the factorizability of the spinon-spinon scattering leads to 

uM = -jUi{e), (64) 

for some constant v. The origin of this relation is recalled in Appendix[B]for completeness. 
There are only two possible values of p at which the scattering amplitudes obey the 
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relation fl64p : p = oo and p = 1. Leaving aside the point p = oo, we conclude that 
integrability at large N requires p = 1. The analysis at this point becomes identical to the 
one carried out by Koberle and Kurak [16] . The large N amplitudes (163|) are then found to 
agree with the large expansion of the minimal refiectionless S-matrix of [H], for which 
n{e) = r2{9) = 0. 

Our analysis above was focused on the scattering amplitudes of spinous, because they 
are the only stable asymptotic excitations for p ^ 1. As far as the integrability is concerned 
it does not really matter whether additional excitations can form in the complementary 
domain p < 1. From a broader perspective it is nonetheless interesting to see what we can 
learn about the spectrum in this domain by considering the scattering amplitudes (1631) . 
We already know that for p < 1 another stable excitation is carried by the gauge field A^. 
It describes a C-odd, SU{N) singlet whose mass is below the two-spinon threshold when 
p < 1. At small enough p this boson is actually the lightest excitation in the spectrum 
since then ~ 12pm^ [3]. It is the large A^ relative of the singlet bound state found 
at A^ = 2, despite the fact that the mass of the latter exhibits a rather different scaling 
~ p^m^ at small p [33]. This feature seems to be tied to peculiarities of the instantons 
gas for N = 2. Other spinon-antispinon bound states will form for p < 1, though, as 
opposed to A^ = 2, none of them becomes degenerate with the gauge field. An analysis 
of the pattern of these bound states based on an effective Schrodinger equation may be 
found in [28J . We can have a glimpse at one of them by mapping the non-relativistic limit 
~ of the scattering amplitudes (!63|) to a Schrodinger equation with a delta function 
potential. One then sees that a shallow bound state is formed for p < 1 in the C-even, 
adjoint channel. Its mass to leading order at large A^ is given by 

^Adj = 2m(l-^^^1^ + ...) . (65) 

We notice that it is at threshold when p = 1. The approximation breaks down when 
p becomes too small. This is not completely surprising since the range of the potential 
increases as ~ small p. Eventually, the interaction becomes confining at p = 0, 

where one expects a rather different scaling for the mass of the bound states [5]. 

Finally let us add a few comments on the large rapidity behavior of the scattering 
amplitudes. In this limit t ~ — s ^ 1 and 6 ~ log(s/m^) ~ Atck/N ^ 1. While U2{0) 
and ri{6) both scale in this case as ~ this is not the case for Ui{6) — 1, which has the 
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following large rapidity asymptotics 




(66) 



This is not exactly what one would expect for an asymptotically free theory and is appar- 
ently related to the fact that we have a circle with finite radius in the UV. Let us try 
to make this connection more precise. It is tempting to believe that the subleading large 
corrections exponentiate in this regime 



For p = 1 this is in agreement with the large rapidity behavior of the finite minimal 
reflect ionless S-matrix, see Appendix [Bl The expression ( 16 7p also coincides with the exact 
S-matrix [30] for N = 2. Thus the equation (1671) seems to be the right guess. If we 
now re-express ( 167|) in terms of the Lagrangian parameter = 2pN we obtain Ui{6) ~ 
exp (— 2i7r/r^). This expression is reminiscent of the large rapidity behavior for the soliton 
scattering phase in the sine-Gordon modelQ Notice that to match the normalization used 
in this paper one should consider sine-Gordon theory for a field ip with period ^/7T. The 
radius would then be given by = /3^/(27r) with being the sine-Gordon coupling. The 
"anomalous" behavior fl^Bl) then seems to indicate that the spinon in this theory is nothing 
else but a soliton for the (p field. This is in line with the discussion in [3] and this analogy 
will appear helpful when considering the theory in finite volume. 

4 Free energy computation 

Having found hints of integrability at large A^, we now wish to verify whether integrability 
is present at finite values of A^. As we pointed out in the introduction, this is of particular 
interest from the viewpoint of the AdS/CFT correspondence, which requires considering 
A^ = 4. According to the analysis of fIE\ and to our previous discussion, the theory at 
p = 1 is described by the minimal reflectionless U{N) S-matrix [T3]. In this section we will 
check whether this assertion holds true at finite A^ and whether there are any modifications 
to the integrability condition p = 1. In order to do this we shall proceed with calculating 

''To make this rigorous one must take into account the proposal of [47] that the sine-Gordon sohton 
scattering phase has to be normahzed as S{d = 0) = — 1. This has the effect of multiplying the Zamolod- 
chikov's S-matrix |12) by a minus sign. 



u^{e) = e-'^/^P + 0(1/6). 



(67) 
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the free energy of the theory at finite chemical potential. This can be done in two different 
ways. One can perform a perturbative computation to the desired order, or make use of the 
conjectured S-matrix. Compatibility of these two computations will allow us, with certain 
degree of confidence, to argue that integrability is present for finite A^. As a by-product of 
our analysis, we will be able to confirm that the value of p, which will be kept arbitrary in 
the perturbative calculation, needs indeed to be set to p = 1 independently of the value of 
A^. 

The computation performed in this section is a standard analysis for two-dimensional 
integrable QFTs. It was introduced in the seminal papers |1HI HH] on the 0{N) sigma 
models and numerous applications to a variety of integrable theories followed, see e.g. 
[Sni EI]- The underlying idea is to consider the ground-state energy density e of a gas 
of spinous with a finite density p. The free energy of interest can then be obtained by a 
Legendre transformation 

f{h) = e{p)-ph, (68) 

with the chemical potential h = de/dp. Computing this quantity directly from the QFT 
leads to an expression in terms of the two RG-invariants of the theory, i.e., the dynamical 
scale A and the deformation parameter p, 

/(/i) = /(/., A, p). (69) 

Thanks to asymptotic freedom, this analysis is tractable at large chemical potential, i.e., 
for h ^ A, where the theory is weakly coupled. On the other hand, the Bethe ansatz 
calculation, based on the asymptotic S-matrix, results in the free energy density directly 
in terms of the mass gap m and arbitrary value of the chemical potential 

f{h) = f{h,m). (70) 

The Bethe ansatz computation will only lead to correct results at those values of p for 
which the S-matrix is a valid physical description of the model. This will be a single value 
for a generic value of A^, while in the case of A^ = 2 the analysis may be done for any p, 
see [33]. The two computations ( 169|) and ( 170|) should match at large chemical potential 
provided that 

m = eA, (71) 
where ^ is a scheme-dependent constant. Since p enters the perturbative computations, it 
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will also get fixed when both computations are compared. 



4.1 Quantum field analysis 

For convenience we perform the QFT analysis using the bosonic formulation (|H]) of the 
f/(l) X SU{N) model. Our starting point is the classical state describing the gas of spinous 
at large density. It is given by the single-spin classical solution 

Zi = e-'^\ zj = 0, j = 2,...,iV, (72) 

where r is the world-sheet time coordinate and u a frequency that we shall soon relate 
to the charge density. This state is interesting because it is directly sensitive to the U{1) 
sector of the theory and hence to the coupling = 47rfi:(l — It has the following energy 
and charge density 

e = k{1 — ri)u'^ , p = 2k{1 — ifjuj . (73) 

We notice that in the CP^~^ limit r/ — ?■ 1, both quantities vanish. This may be easily 
understood. In the CP^~^ model the state ( 172|) is equivalent to the Goldstone vacuum up 
to a gauge transformation and has therefore zero energy. For the free energy density we 
have 2 2 

f{h) = -K{l-r^)h' = --^, (74) 

with the chemical potential h = de/dp = u. 

Let us proceed now to the one-loop analysis. It is convenient to go to the Euclidean 
space r — )■ —icro- To compute quantum corrections to the free energy, we shall expand 
the Lagrangian around the solution (172|) and calculate the partition function of the theory 
Z{h), which will allow us to extract the free energy 

A convenient parameterization around the background solution fl72|) is given by 

z, = e-'^^==^, ^, = 6-"'^'^-=^, J =2,..., AT, (76) 

where ha = h^a^, with = {h,0) and = ((To7<^i) the Euclidean world-sheet coordi- 
nates. This parametrization breaks the symmetry down to U {N — l)xU{l). Note that we 
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used that 00 = h to have the chemical potential as a "background field" explicitely. More 
importantly, we have also allowed for a chemical potential ah for the — 1 coordinates 
Zjj^i. This does not change the semiclassical limit since perturbatively Xj ~ Vv^ ~ 0- 
introduces however a new parameter a, which eventually must be eliminated. This will be 
done by demanding that the system is in its ground state at fixed h. To see how it can be 
done, we plug the expressions (176|) into the Lagrangian and expand up to quadratic order. 
We find up to total derivatives 

£ = -,^(1 - r/)/i2 + £2 + . . . , (77) 

with 

C2 = K {d^ + I3h^) X {d^ - X + <1- vfh^XX + '^(l - V)d,^d^^ ■ (78) 

The new parameter (3 = a — rj. We see that the fields Xj have acquired a mass (1 — ri)h, 
while the field has remained massless. The fields Xj further couple to the chemical 
potential For /3 7^ the system therefore forms a condensate of massive bosons at 

rest. Minimizing the energy then amounts to choosing /3 = 0, or a = 77. The fields Xj 
become free and their contribution to the free energy density is easily evaluated. The 
massless field does not contribute if one uses dimensional regularization. 

The evaluation of the one-loop determinant for the — 1 complex bosons is straight- 
forward. It leads to 

/(^)-^h^('-(^;:i)4)-H' 

when expressed in terms of the radii. Along the way we absorbed the UV divergences into 
the renormalization of the radius r, 

+ ^^^^^ (^-^^ - log47r + 7e) ^ , (80) 

here performed in the MS scheme and with 7e the Euler-Mascheroni constant. This re- 
definition is in agreement with the renormalization group equations iQ derived from |29] . 
This is most easily seen at the level of the free energy density (1791) . The latter defines a 
physical observable and its dependence on fi should drop out when the couplings R^,r^ 
fulfill the RG equations This is easily verified to be the case, at the given order in 
perturbation theory. We can also check the consistency of the expression (179|) at special 
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points. In the 0{2N) limit, when = R"^, we obtain 



fo{2N){h) 



1 + 



2(A^- 1) 

^2 



log 



^1) 



This result agrees with the one of [IH] when the couplings are appropriately related. In 
the opposite limit, — )■ 0, we find that f{h) vanishes, as expected for the CP^^^ model. 

To compare with the S-matrix computation one needs to express the free energy density 
f l79p in terms of the RG- invariant parameters of the model, i.e., in terms of the dynamical 
scale A and deformation parameter p. Using the expressions f lTUj) and f ITT]) we arrive at 



f{h) 



ph'^N 
2tx 



1 - 



p{N - 1) 
A^log(VA) 



p{N -1){N + P-2) log log (/i/A) + pN{N - l)d 
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m log'(VA) 
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now understood as an expansion valid at large chemical potential h ^ K. The constant 
djq reads 



cLn = - \ogp + 



1 



2p 3p-2 
~ ^ 2N 



(83) 



4.2 Matching the S-matrix analysis 

To calculate the expression for the free energy of gas of spinous using the conjectured 
S-matrix, we first observe that only the spinon-spinon scattering phase in the symmetric 
channel, S{9) = Ui{9) + U2{9), is relevant for that purpose. The explicit expression for S{9) 
may be found in Appendix [Bl It turns out that S{9) is identical to the scattering phase for 
fundamental excitations in the SU{N) chiral Gross-Neveu model, up to the substitution 
1/A^— j-l — 1/A^. The computation of the free energy density in the SU{N) chiral Gross- 
Neveu model was carried out in [50j. One can thus directly translate the result to the case 
at hand. For h ^ m we find 



N 



{N-lf\og\og{h/m) + N{N-l)D 



N 



N log (h/m) 



m \og^h/m) 



+ ... 



where m is the mass of a spinon. The constant D^^ stands for 

log2 1 3 



Dn = \ogT{l + l/N) 
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2 ^ 2A^ 
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We immediately notice that the logarithmic pattern of the QFT and Bethe ansatz 
expressions, Eq. fl79|) and f lH^ respectively, are compatible. They match precisely if and 
only if p = 1 for any A^. It confirms that integrability occurs at p = 1 without any 
corrections. One other inference that may drawn is that the exact S-matrix has 
been properly identified, else the structures of both expansions would be different. Fixing 
further terms allows us to relate the mass gap m to the dynamical scale A, i.e., determine 
the constant ^ in (17T]) . It is given by 

m_ (2/e)V^ 
A -r(l + l/Ar)- 

Note that this formula is valid in the MS scheme. At large A^ we immediately verify that 
m = A + 0(1/A^), as required by the large A^ analysis. For N = 2, the relation between m 
and A is known for any p 



"«=2(rt = (2/e)'-/^ei±iV!)A. (87) 

Tip 

It is easily seen to agree with (!86|) when p = 1. We note also that for p = 2 the mass 
gap f l87|) is identical to the one of the SUSY CP^ model [ST], in line with the discussion in 
Sectional 

Finally, we observe that the expression (1861) is in agreement with the result reported 
in [28]. This reference did not include any details on its derivation however. Moreover, 
it apparently overlooked the renormalization of the radius r, as pointed out in [29J. The 
above analysis shows that the running of this coupling is actually necessary for a proper 
match of the QFT and S-matrix computations. We notice nevertheless that this feature 
is specific to our choice of probe. The computation of the constant ^ could also be done 
by performing the one-loop analysis around a different classical solution, for instance the 
two-spin solution ~ z(zj^. The free energy corresponding to the latter solution would not 
be sensitive to the renormalization of the radius r, which could therefore be treated as a 
constant. 



5 The AdS/CFT effective model 

In this section we analyze in more details the particular case of the CP^ model coupled to 
a Dirac fermion with charge k = 2, which was proposed in [20] to govern the low-energy 



33 



effective theory of the Gubser-Klebanov-Polyakov (GKP) string [22] in the long-string hmit. 
For reader's convenience we recall the form of the Lagrangian 



£ = K{d^ - iA,)z{d^' + tA^^)z + i^Y{d^ - 22A^)V^ + ^ (V^7mV^)' • (88) 

According to our previous analysis the fermonic model ([2]) is integrable when the UV value 
of the Thirring coupling fulfills 

Aoo = 7r(A;VAr - 1) = , (89) 

where in the last equality we have used the parameters of the string theory model. Please 
note that the absence of the UV value of A goes in line with the form of fISS]) . The running 
of the Thirring coupling was found in fl29|) to be given by 

which upon plugging back into ([2]) leads immediately to (IHHj) ! It is satisfying to find that 
the Lagrangian f l88|l . derived from a renormalizable and UV complete superstring action, 
displays the correct one- loop induced value of A, which is scheme independent. We may 
thus conclude that ( 155]) is integrable and belongs to the class of integrable U{1) x SU{N) 
models. This corroborates the conjecture that the original AdS^ x CP^ string theory sigma 
model is quantum integrable. Since this sigma model has a gauge theory dual and a set of 
all-loop spectral equations was conjectured for the said duality [25], we will independently 
propose Bethe equations for fl88|) . 



5.1 Physics in finite volume 

In this subsection we consider the fermionic model on a cylinder of length L. We want 
to address the problem of constructing the spectrum in the asymptotic domain mL ^ 1, 
where m is the mass of a spinon. To do that we will need to elucidate the role played 
by the fermion charge k. We will assume it to be integer, while the value of the Thirring 
coupling will be fixed requiring the integrability condition to be fulfilled. The spectral 
problem relevant to the string theory will then appear as a special case of = 4 and 
k = 2, and shall be discussed in more detail later on. 

It is well known that for integrable theories the large volume spectrum is encoded in 
a set of asymptotic Bethe ansatz (ABA) equations. There are essentially two pieces of 
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information required to write down such equations. The first ingredient is the S-matrix, 
which for the model at hand may be found in Appendix [Bl The second one is the choice of 
boundary conditions for the multi-spinon wave function. To minimize the amount of tech- 
nical details, we will first look at scattering of K spinous of the same SU{N) polarization. 
The S-matrix in this case reduces to the phase S{6) introduced before. The corresponding 
ABA equations, which are the quantization conditions for the momenta of the spinous, 
take the following form 

K 

^-^V{0,)L ^ ^ JJ g^Q^ k = l,...,K, (91) 

where p{6) = msinh6' is the momentum of a spinon with rapidity 6. The phase q is 
associated to the monodromy of the wave function as one of the spinous goes around 
the cylinder. Assuming that the bosonic fields of the theory are subject to periodic 
boundary conditions, one would naively choose q = 1- We shall argue below that this is 
not the proper choice for generic values of k. 

So far the fermion charge k was subsumed into the parameter p and did not play a role 
on its own. We have, however, already mentioned its relation to a discrete = 'Z2k/'^2 
symmetry [5]. Not surprisingly, understanding the effect of the charge k on the ABA 
equations parallels the implementation of this symmetry. To clarify this point we shall 
first develop a useful analogy with the fc-folded sine-Gordon (SG) model 152]. 

The latter model can be defined by the Lagrangian 

^SG = ^d,^d^^ + ^ COS {2kV^^) , (92) 

where ^9 is a compact boson with period a/tt. Forgetting for a while the SU {N) symmetry of 
our problem, the model (l92i) can be seen to have a lot in common with the model studied 
in this paper. They share for instance the same symmetry which is generated by 
+ ^fii jk. Moreover, in both cases this symmetry is spontaneously broken in infinite 
volume. In the SG theory this is directly observable at the level of the Lagrangian (l92i) . 
The theory has k degenerate vacua in a given period, which are located at <y9 = Uy/Ti / k with 
?T, = 0,...,A; — 1. They are all equivalent and each breaks the discrete symmetry. In our 
case, the phenomenon is not visible at the tree level nor at any finite order in perturbation 
theory. It can be revealed however following the observation made in [3J that spinous are 
solitons from the perspective of the bosonized fermion ip. In the background of a spinon 
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the field jumps by 



(93) 



where Q = 1 is the spinon U{1) charge. This immediately follows from the Gauss law ( 126|) 
upon bosonization, or equivalently from (IT^ . Since the spinon is a stable massive particle 
carrying the minimal amount of U{1) charge, one expects degenerate vacua separated 
from one another by Aip = ^fiilk. Up to inessential details like the explicit form of the 
potential or the expression for the soliton S-matrix, this picture is identical to the one 
emerging from fl^?2]) . 

The charge k is associated to rescaling of the field <y9 and may always be eliminated in 
infinite volume. In finite volume, on the other hand, the function of the charge k becomes 
somewhat more manifest. For definiteness we will assume periodic boundary conditions 
for the field i.e., ^{a + L) = v'(cr) mod a/tt. Due to tunneling effects the k vacua 
are no longer localized and their energy degeneracy is lifted leading to restoration of the 
Zfc symmetry. Moreover, the entire Hilbert space of the theory splits into k subsectors 
associated to different representations of the symmetry. Each subsector is associated to 
a Bloch wave with quasi-momentum P„ = 2Tm/k. The effect of a Bloch-wave background 
on the ABA equations is known [22] '■ the wave function acquires an additional phase 
shift g„ = exp {iPn) each time a soliton circles around the cylinder. An anti-soliton picks up 
the inverse phase l/q-n- One way of understanding it is by recalling |13] that a Bloch-wave 
background is equivalent to inserting the vertex operator 



as required for a Bloch wave with momentum P„. This background acts non-trivially on 
excitations with non-zero topological charge or winding number W = A(p/-\/7r such that 
they pick up a phase exp {ikWPn) when transported once around Vn- This applies in 
particular to a soliton. The latter, interpolating between two adjacent vacua, contributes 
W = 1/k resulting in the phase shift g„. 

The overall effect of the charge k on the ABA equations is thus the division into 




(94) 



at the bottom of the cylinder. Under the shift — + ^pn jk it transforms like 



(95) 
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subsectors characterized by a twist g, which itself is a root of unity 




(96) 



There is a further modification \5T\ . The number of sohtons minus antisohtons is quantized 
in the units of k. This selection rule originates from the quantization of the total winding 
number = (v5(cr) — iflcx + L)) / ^/^T . Since a soliton contributes only 1/k unit to this number, 
at least k solitons are needed to make it integer. 

The above remarks clarify why the charge k is related to the ABA equations. It is 
not yet clear however how exactly it modifies the spectrum of the fermionic theory. We 
shall see that equation f ^B]) has to be slightly corrected. In fact equation is pertinent 
to the Zfc quotient of the bosonic model. This claim might seem puzzling given that we 
have argued before that the bosonic and fermionic models are equivalent. Nevertheless, 
examples of theories are known, which, despite being equivalent on the infinite plane, have 
different spectra in the finite volume. The most celebrated example is given by the duality 
between the sine-Gordon and the massive Thirring models |17]. The difference between the 
two sets of Bethe equations is tiny but significantly changes the spectra in finite volume. 
Here we are facing a similar problem. 

There is a matter in which the bosonic and fermionic theories differ from each other. It 
is the action of the Z2 symmetry (f ^ (f + a/tt, as mentioned in Section 12.31 This operation 
is trivial in the bosonic theory, where the field ip is A/vr-periodic from the beginning. In the 
fermionic theory, however, this transformation is associated to the operator (—1)-^, with 
the fermion number 



The corresponding spaces of local (gauge-invariant) operators of the two theories are dif- 
ferent due to the presence of operators anticommuting with (—1)^ in the operator space of 
the fermionic theory. This is the root of the difference in the ABA equations for the two 
models. We refer the reader to [17] for a more detailed discussion in the context of the 
sine-Gordon / massive-Thirring duality. 

We can now derive the correct twist q for the spectral equations of the fermionic theory. 
We assume antiperiodic boundary conditions for the fermion ip while the bosons are taken 
to be periodic. There are two reasons for this particular choice. Firstly, we believe it is the 
proper set of boundary conditions to describe the string theory spectrum corresponding to 



F = number of fermions — anti-fermions . 



(97) 
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= 4 and k = 2. Secondly, the Neveu-Schwartz (NS) sector is the simplest one from the 
perspective of the state-operator mapping, which is a one-to-one correspondence between 
finite- volume eigenstates in the NS sector and local gauge- invariant vertex operators of the 
theory. These vertex operators are typically of the type 

V ~ VfVb , (98) 

where Vp and Vb are written in terms of fermionic and bosonic fields, respectively. The 
bosonic part has a more transparent interpretation in terms of spinous. For example, the 
operator Vb ~ z*^ . . . z*-^ is likely to correspond to a finite-volume state made out of K 
spinous in a totally symmetric representation of SU{N). What needs to be understood is 
the effect of the fermionic component Vp on these spinous. The fermionic operator has an 
electric charge which blots out the U{1) charge of the bosonic operator Vb- This leads to 
a selection rule on the total number of spinous. We will come back later to this issue. The 
important point is that the vertex Vp has also a "Bloch-wave" in close similarity to the 
operator Vn considered in flMl) . According to the rules of the bosonization, the operator 
Vp with fermion number F should contain the wave 

exp (^iriy^ip) , (99) 

where n is an integer fulfilling (—1)" = (—1)'^. By analogy with our previous discus- 
sion, we expect therefore that a spinon with winding number 1/k will pick up an extra 
phase exp (inn/k) due to the fermionic "background". This heuristic argument leads us to 
propose 

q' = {-ir, (100) 

as the phase for the fermionic theory in the NS sector. We observe that for F even it 
reduces to the twist f l96p for the bosonic model. In other words, the bosonic and fermionic 
ABA equations are the same for states which are neutral under the Z2 symmetry. This is 
analogous to what happens for the sine-Gordon / massive-Thirring duality [17j, which in 
our notation corresponds to k = 1. We also notice that fllOOj) no longer defines a twist in 
the original sense since q depends via F on the state considered. Finally, we should impose 
a selection rule for the fermionic theory. It may be traced back to the gauge invariance 
of the theory, which requires that the total bosonic U{1) charge Q = K — K, with K the 
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total number of anti-spinons, is a multiple of the fermion charge k. Explicitly, 

K~K = kF. (101) 
This selection rule is the same as for the bosonic theory, since in the NS sector F G Z. 

5.2 ABA equations 

We are now in the position to present the complete set of ABA equations for the NS sector 
of the string model (155]) . To resolve the mixing related to the SU{A) symmetry, one has to 
diagonalize the monodromy matrix associated to the minimal refiectionless S-matrix. The 
analysis can be done by means of the algebraic Bethe ansatz and is very similar to the one 
performed for the alternating Heisenberg spin chain of the planar ABJM theory We 
refer the reader to this reference and here quote only the final result 

j=i j=i I 

K r\n / \ i ^1 ■ ^2 i 

nui^k — ^C7j/7r + 2 _ TT ""1,*: ~ ""iJ + * TT ^^''^ ~ ~ 2 
, -oft.u - i. ~ 11 „, _ 7,, , _ 11 ' 



K2 . ■ Ki i Ki i 

}4 - «2,i - « U2,k - U3J + \ y U2.k - Mlj + f ' 



, _ 9fl. _ i ~ 11 , _ , _ 11 



i ' 



i/q n ^(^"^ - ^".) n ^1 (^"^ - n '''' ^ 



jjLk j=l 



f \ 2ek/TT - na,, - f 



In these equations S{6) and ^1(6*) are the spinon-spinon and spinon-antispinon scattering 
phases, see Appendix [B] for further details. The momentum of a spinon or anti-spinon, 
carrying respectively rapidity 6 or 6, is given by p{6) = msinh6' and similarly for the 
anti-spinon. The total energy is given by 



K K 



E = m cosh 9k + m cosh 9k ■ (103) 
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The numbers count the isotopic roots Mi,M2 and M3, which change the flavors of the 

K spinous and K anti-spinons, as may be inferred by looking at the su(4) Dynkin labels 
of the state 

gs] = [K- 2K^ + K^, K, + K, - 2K2, K - 2K, + K^] . (104) 
The twist q is solution of 

g' = (-ir, (105) 
which, together with (llOip . allows one to write 

g = ± exp {i7r{K - K)/A) . (106) 

We stress that both solutions have to be considered. In the F-odd sector this is required by 
the parity invariance of the theory, i.e., if (6*^, ^fc, Uj,fc) is a solution of the ABA equations 
then so is {~0k, —Ok, —Ui^k)- This property is not guaranteed if g 7^ l/q. Finally, we 
notice that an immediate consequence of the ABA equations is that the total momentum 
P = Y^kPi^k) + Y^kPi^k) satisfies 

e-^P^ = q'<-K = {-If . (107) 
This is in agreement with our choice of boundary conditions for the string model. 



5.3 State-operator matching 

This subsection is devoted to investigating in more detail the mapping between solutions 
to the ABA equations (I102p and local operators of the theory. According to this correspon- 
dence a state of energy E and momentum P is associated to a vertex operator of scaling 
dimension A and spin S. This relation should become more and more evident at small 
length mL <^ 1, for which IMl ES] 

^ TTC 27rA , , , 2713 

E = -— + ^ + o{l/L), P = -^, (108) 

with the UV central charge c = 2 —1 = 7. The second equality is actually expected to be 
valid for any L. The difficulty that one immediately encounters at the attempt of verifying 
(11081) with the help of fll02p is that the ABA equations are only approximate. They do not 
take into account off-shell processes like vacuum tunneling, exchange of virtual particles. 
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etc., which become more and more important at smaller values of L. As a consequence, 



the ABA equations often fail to reproduce fllOSp making the mapping between operators 
and states hard to quantify. Despite these complications we will present some evidence 
that the twisted ABA equations correctly capture the properties of the low-lying energy 
eigenstates. 

Let us start with the two "vacua". These states are easily found: they correspond 
to no root at all and have exactly zero energy in the ABA description. This is a trivial 
observation but it is nevertheless in agreement with the expectation that at large volume 
the vacua are exactly degenerate. The degeneracy should be lifted by tunneling processes 
as in the folded sine-Gordon theory [121 This is confirmed by the Liischer formula, 
which corrects the ABA equations at large volume mL ^ 1. With the twist q included 
the finite-size correction writes 

E„ ~ -—Ki{mL) - —Ki{mL) . (109) 
TT Tiq 

Here, Ki{z) is the modified Bessel's function with the asymptotics Ki{z) ~ for z ^ 1. 
We derived f llOQp by adapting to our case the analysis performed for the Bloch wave vacua 
of the sine-Gordon theory |13] such that the = 4 flavors of spinous and anti-spinons are 
correctly incorporated. The two NS vacua correspond to the two possible choices g = ±1. 
The formula f llOQp shows that the degeneracy is lifted and that the true vacuum has q = 1- 
The NS ground state is thus Z2 even, as expected. It corresponds to the identity in the 
operator picture. But what is the vertex operator that creates the uplifted vacuum? It 
has to be odd under the Z2 symmetry — > + \fTi jl and singlet under S'f/(4). Moreover, 
among all such operators, it is likely to be the one with minimal scaling dimension A in 
order to minimize the energy gap at small volume 

E_-£;+ = ^ + ... . (110) 

There exist two possible candidates, 

~ cos (2-\/vr(y9) , and '^^•o^ ~ sin (21/77^9) . (HI) 



They are both spinless and have dimension A = 1 to leading order at weak coupling 



''Perturbative corrections to the scaling dimension are irrelevant here. Thanks to asymptotic freedom 
they are associated to subleading contributions at small length. 
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To decide between the two we will look at the bosonized forms of the operators and use 
some basic quantum mechanics. We see that only the first operator describes a state 
with a wave function peaked around the (expected) minima = 0, A/vr/2 of the potential. 
On the contrary, the second operator depicts a wave function localized at the tips of the 
potential and is more likely to be an excited state in infinite volume. In any case, one of 
the two operators, or a linear combination perhaps, has to be an excited state and should 
be embedded as a non-trivial solution in the q = —1 subsector of the ABA equations. 
Assuming a singlet of SU{4), the simplest solution in this subsector has vanishing roots 
6 = 6 = ui = U2 = U3 = 0. Interestingly, its energy in the ABA approximation is exactly 
2m for any length L leading to a suggestive interpetation: it may be seen as the C-odd 
singlet bound state at rest which we know is exactly at threshold. This interpretation is 
in line with the comments in [17] on the vertex operators for breathers in the sine-Gordon 
theory. 

A more direct test of the state-operator correspondence is available if one looks into the 
non-singlet U{1) sector. Due to the selection rule the simplest configurations correspond to 
{K, K) = (2, 0) and its charge conjugate. Their fermion numbers are F = ±1, respectively, 
and the allowed values of the twist are g = ±z in both cases. The solutions to the ABA 
equations are easily constructed at small length since in this regime the rapidities are large 
and we only need the asymptotic expression for the scattering phase S{9), see Appendix [Bl 
The solutions with F = 1 and q = ±i, which minimize the energy, are given to leading 
order by 

Pi-±^ P^-^±^- (112) 
The total energy and momentum take the following values 

They should be compared with fllOSp after subtracting the energy of the vacuum (~ 
— 7rc/(6L)) and using the labels of the associated vertex operators. These should have 
fermion number F = 1, belong to the symmetric representation [2, 0, 0] of S'f/(4), and have 
minimal scaling dimension. There are precisely two gauge-invariant operators with these 
properties, 

ip±z'z^ , (114) 

where tp± are the spin S = ±1/2 components of the Dirac field. Their scaling dimension is 
A = 1/2 to leading order at weak coupling, in perfect agreement with ( 11131) and fllOSp . It 
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is amusing to note that the energy and momentum of the states flll3p come directly from 
the dimension and the spin of the Fermi field in the operator picture f lll4p while they are 
distributed over the two spinous in the ABA description flll2p . Similar comparisons can be 
performed by considering more excited operators of the type ilj±D±ip±z^z^ z^z'' , etc., with 
D± being the light-cone covariant derivatives. 

This concludes our tests of the twisted ABA equations fll02p . Our discussion provides 
a somewhat appealing, though not very strong, evidence for the correctness of these equa- 
tions. To put the analysis on a firmer ground, it would be interesting to write down and 
analyze the full-fledged thermodynamical Bethe ansatz (TBA) equations for the string 
model in the NS sector. As opposed to the ABA equations, the TBA equations should be 
valid for any length L. An analysis of this type was performed for the fc-folded sine-Gordon 
theory in [52]. 

6 Conclusions 

Integrability is a rare and unique property of quantum field theories. It is usually difficult to 
prove it rigorously since this would require a good handle on the non-perturbative physics 
of the theory. There are however certain quantities and features, one can study using 
perturbation theory or some non-perturbative methods, that provide hints as to whether 
the quantum integrability is present. 

In this paper we have studied a family of U{1) x SU{N) theories, which may be formu- 
lated either exclusively in terms of bosonic degrees of freedom or by coupling a massless 
fermion with non-zero self-interaction. Both formulations are equivalent quantum me- 
chanically and are equipped with a continuous parameter p. It is related to the UV values 
of the coupling constants of the bosonic and fermionic models by f fT2|) and f l25|l . respec- 
tively. Since by definition it is RG-invariant, it can be deemed physical at the quantum 
level. The corresponding {p, N) family of models covers a wide range of theories, which 
are summarized in Table [H 

We have found several indications that the class of models considered may be rendered 
integrable with the right choice of the physical parameter. This fine-tuning is not necessary 
for = 2 because the U{1) x SU{2) model is equivalent to an integrable one-parameter 
deformation of the PCF model. Our analysis suggests that the theory should remain 
integrable at higher values of at least for p = 1. This may be inferred already at the 
level of the large scattering matrix, where integrability places severe constraints on 
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p = 


< p < oo 


p = CO 


N = 2 


0(3) model 
+ free boson 


U{1) X SU{2) model; 
integrability at any p 


0(4) model 


N>2 


+ free boson 


U{1) X SU{N) model; 
integrability at p = 1 


0{2N) model 



Table 1: Different limiting cases of the model under consideration. 



its structure. A similar conclusion should be reached by studying the 2 — > 4 amplitude 
at leading order in large N. An important question is whether this choice of p receives 
any modification beyond the large N limit. The one-loop free energy computation in 
Section m is a direct confirmation that this is not the case. This leads us to believe that the 
f/(l) X SU{N) model is integrable for p = 1. We should stress however that this might not 
be the only point where integrability prevails. Alternative finite values of p might exist, 
though presumably only at small N, and it would be of interest to investigate more closely 
this eventuality. 

A model belonging to the integrable class identified in this paper has been recently 
found in [20j . It governs the dynamics of massless excitations around the GKP solution of 
the ^4^5*4 X CP^ string theory. These excitations comprise the CP^ degrees of freedom and 
a Dirac fermion. The Lagrangian of the model is embedded in the fermionic formulation 
(|2]) and its parameters fulfill the integrability condition p = 1. The integrability of this 
effective model is a strong evidence in favor of the integrability of the full AdS^ x CP^ 
super sigma model. In Section O we put forward a set of asymptotic Bethe equations for 
the fermionic model. Our argumentation relied solely on the physics of the model and was 
not inspired by the AdS/CFT correspondence. This set of spectral equations may thus 
be used to test the veracity of the all-loop ABA equations proposed for the AdS4/CFT3 
duality 
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A Renormalization of the Thirring coupling 



In this appendix we show that the Thirring couphng runs. In order to do this, we will 
construct the two-point function for the fermionic current = ip'^^ip up to order 0(1/k^) 
in the weak coupling expansion. More precisely, we shall consider the form factor D{p^) 
defined by 

d^x {nx)fm e^'^ = -^^^ (^V'- - ^) . (115) 

Notice that this quantity is dimensionless in two dimensions {D = 2). An important 
property of this correlation function is that it is observable and hence UV finite. We 
will demonstrate that the finiteness cannot be preserved if the Thirring coupling does not 
renormalize. 

We should first comment on the renormalization of the bosonic coupling l/n itself. At 
leading order the running of this coupling constant follows that of the CP^~^ model. The 
reason is that the contribution from a fermionic loop are controlled by the correlator flllSp 
evaluated in the Thirring model. The latter is well-known to be finite, 

^(p')Thirnng = 1 + A/vT + 0{D - 2) , (116) 

for D = 2. Hence all divergent corrections to the coupling 1/k, come from bosonic loops 
at leading order. In the MS scheme the renormalized coupling 1 / Kj-e^ is then defined as in 
the CP^-i model 



N 

4:71 



2 

7e - log 47r 



D-2 



+ 0{1/k). (117) 



We refer the reader to [H [3] for an explicit derivation of this relation and for a more detailed 
discussion of the renormalization of the CP^~^ model. 

In order to calculate D{p'^), we shall first derive the effective propagator for the gauge 
field A^, 

^ dPx e^P^ l^N'{x)A\^))cp^-^ = D^'^{p^) , (118) 



by integrating out the bosonic degrees of freedom of the CP^^^ model. To leading order 
at weak coupling, this model is described by free massless complex bosons coupled to a 
gauge field. The relevant part of the Lagrangian is 

C = Kd^xd^X - ^^^A'^x'^X + ^^A^A^ + . . . , (119) 
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Figure 4: Gauge field propagator to lowest order. 

+ i^^-"-*--^ 1^ + • • 



Figure 5: Geometric series for the diagrams contributing to the divergent part of the 
current-current correlator. Crosses stand for insertions of the fermionic current. Dashed 
lines indicate contractions performed in the massless Thirring theory. 

^ y 

where x = iX2, ■ ■ ■ ,Xn) is a multiplet of — 1 complex bosons and xdfj,x = xd^x ~ 
d^XX- Note that these fields are not subject to any constraint since they parameterize the 
transverse directions around the Goldstone vacuum z = (1, 0, . . . , 0). The contribution of 
these massless bosons to the propagator of the gauge field is depicted in Fig. HI A direct 
computation of the one-loop diagram leads to 



^ 2k 4^2 ^ J (27r)^ (p + g)V V / V ; 

where the scale fi was introduced to keep the coupling k dimensionless. A straightforward 
algebra yields 

2„z> f {p + 2qY{p + 2qr Ijp') f p^p^\ 



(27r)^ (p + g)V Z^-lV P 



where 



' y (27r)^(p + g)V 2^sin(7rZ}/2)r(D/2-l/2) V ' 

Notice that this integral is divergent for D = 2. 

We are now in the position to compute the correction to the two-point function f lllSp . 
or more precisely the divergent part thereof. We calculate it up to ©(I/k^) in the weak 
coupling expansion, but to all orders in the Thirring coupling. The relevant diagrams are 
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Figure 6: Example of a diagram contributing to the finite part of the correlator. It rep- 
resents the connected part of the four-point function of fermionic currents in the massless 
Thirring theory convoluted with the gauge field propagator. Since the first one is of order 
0{D — 2), it compensates the pole in the latter ~ l/{D — 2) leading to a finite contribution. 

shown in Figure O Their contribution is a geometric series that sums up to 

D{p' = ! + - + - / V + ^^/'^ • 

We stress that this computation may not be accurate enough to capture the finite part 
of the perturbative correction ~ The reason is that we are disregarding certain 

diagrams which are superficially of order 0{D — 2) but might nonetheless make a finite 
contribution when convoluted with the divergent propagator (I120p . An example of such 
diagrams is depicted in Figure |6l Our analysis should reproduce correctly the pole at 
D = 2 however. Taking the limit D ^ 2 \\i f lT23|) with help of flT22|) we get 

D(p^) = IH \ -) f- + ... . (124) 

The result is UV divergent. The part linear in N will be made finite after renormalization 
of the bosonic coupling, see flll7p . To absorb the remaining divergence we must renormalize 
the Thirring coupling, 



A+ — 



7e - log Atx 



^A, (125) 



LP -2) 

where we again work in the MS scheme. The resulting RG equation takes the form 

This is the expression we used in Section 12.31 It was found to be in agreement with the RG 
equation of the dual bosonic theory. The calculation performed here is not very different 
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from the free energy density computation done for the bosonic model in Section HI Both 
lead to the same conclusion that the deformation parameter of the bosonic and fermionic 
model must be renormalized. 



B Exact S-matrix 

In this appendix we discuss some properties of the factorized S'-matrices with U (N) sym- 
metry. Their complete classification has been proposed long time ago by Berg et al. [H]. 
It was found that there are several classes of solutions and the minimal representatives of 
each class do not depend on continuous parameters for > 1. This immediately implies 
that there is no continuous deformation of the minimal 0{2N) S-matrix [12] preserving 
the f/(l) X SU{N) symmetry^ This important result can be used to argue that the inte- 
grability of the U{1) x SU{N) models considered in this paper can only be sporadic with 
respect to the deformation parameter p. 

The absence of a continuous deformation of the 0{2N) S-matrix is however in disagree- 
ment with the S-matrix derived by Wiegmann [30! for N = 2. He found a solution given 
by a tensor product of the minimal SU{2) S-matrix and the sine-Gordon S-matrix. The 
latter contains a free continuous parameter p. In this appendix we would like to clarify 
this apparent paradox. Contrarily to what has been stated in [T3|, we show that the case 
= 2 is special and permits the solution found by Wiegmann. The space of solutions for 
= 2 is enlarged in the very same way as in the case of 0{N) factorized S-matrices [T2] . 
Explicitly, we observe a reduction of the number of equations constraining the S-matrix 
due to additional identities inherent to = 2. We have no doubt that this caveat has 
been known before. We did not find however any discussion nor comment in the literature, 
which is why we decided to comment on it below. 



B.l The Yang-Baxter equations 

The equations that implement the factorization of the S-matrix are the so-called Yang- 
Baxter equations. Their construction is facilitated by the use of the Faddeev-Zamolodchikov 
algebra, see for instance [12]. The spinon (antispinon) carrying rapidity 6 is denoted by 
Zi{6) (respectively Zi{6)) with the SU{N) index i = 1,...,N. The incoming states are 
written as strings of such symbols ordered with decreasing rapidities. An opposite ordering 

^The 0{2N) S-matrix was found to be embedded in the class III of the classification, see discussion 
in [Hj. 
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is used for the outgoing states. Their algebraic relations are encoded into the commutation 
relations between different Z's and/or Z's. For example the spinon-spinon scattering is 
implemented by 

z,{e,)z^{e2) = s^j{9)Zii92)Zkie^) , (127) 

where = 6i — 62,61 ^ 62, and S^j{6) is the two-to-two spinon S-matrix. Similarly, for the 
spinon-antispinon scattering one has 

Z,{6i)Z,{62) = Fg{6)Zii62)Zk{6i) + BI'; (6) Zj, (62)21(61) . (128) 

The matrices F and B stand for the forward (transmitted) and backward (reflected) compo- 
nents of the scattering matrix, respectively. The corresponding formulas for Z are obtained 
by charge conjugation. 

Imposing the U(N) symmetry allows to decompose the matrices S, F, B into six scalar 
amplitudes, as done in fl58|) . The crossing symmetry yields three linear relations between 
them, see (l59i) . Unitarity, on the other hand, provides quadratic constraints. These rela- 
tions are most easily written in terms of the amplitudes associated to the various SU (N) 
invariant channels [H]. By way of illustration, after introducing the scalar factor for the 
scattering of two spinous in the symmetric channel, 

S(6) = Sll(6)=m(6)+U2(6), (129) 

one arrives at the unitarity equation S(6)S(—6) = 1. There are six such equations for the 
six 5'f/(A^)-invariant channels of spinon-spinon and spinon-antispinon scattering. Their 
explicit expression can be found in [14j. 

Finally, we have the Yang-Baxter equations. These are cubic in the S, F, B matrices 
and take care of consistency of the factorization of the three-body scattering [12]. They 
are equivalent to the associativity of the Faddeev-Zamolodchikov algebra. In our case we 
need to consider the four distinct processes with initial and final states listed below 

Z,(6i)Z,(62)Zk(6,) ^ Zn(6s)ZU02)Zi(6i) , 
Z,(6i)Z,(62)Zk(63) ^ Z„,(6,)ZUe2)Zi(6i) , 
Z,(6,)Z,(62)Zk(6,) ^ Zr,(6,)ZU02)Zi(6,) , 
Z,(6i)Z,(62)Zk(6s) -> Z„,(6,)ZU02)Zi(6i) . 

The inequivalent ways of computing the S-matrix for these processes, obtained by iterative 
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use of the commutation relations, yield cubic matrix identities between S,F,B. The equa- 
tions corresponding to the first two lines in (11301) constrain the spinon-spinon scattering 
matrix S and relate F to S. They are oblivious to the backward scattering matrix B. Each 
of them provides a single functional relation, 

fiOl2) + /(^23) = /(^IS) , h{e,s) - hie^s) = -/(^is) , (131) 

where Oij = 9^ - 0jJi0) = Ui{e)/u2i9), and h{e) = ti{e)/t2{e). The general solution 
to these relations is well-known [H] . After taking into account the crossing relations 
between ui^2{d) and ^1,2(6') it reads 

ui{e) = '-^U2{e) , tr{e) = -'-{o- m) t2{9) , (132) 

where the parameter u is still to be determined. The relations fll32p apply to all U{1) x 
SU (N) factorized S-matrices which are of interest to us, that is assuming that Ui{6),ti{6) ^ 
0. 

For theories with reflectionless scattering, when B = 0, the relations f ll32p are all we 
need. In all other cases one should also analyze the equations for the last two processes in 
f ll30p . They relate the backward scattering B to S and F. We present these equations in 
a pictorial form in Figure [3 



B.2 Minimal reflectionless S-matrix 

We will first turn to the minimal reflectionless S-matrix since it is of particular significance 
for the analysis carried out in this paper. By definition it is characterized by ri, r2 = 0. In 
this case the unitarity equations fix the value of the parameter u in f ll32p to 

and lead to ^ 

ti{e)h{-e) = 1 , u{t7T - oMzTT + e) = -J—^ . (i34) 

t/"^ + 

The minimal S-matrix is a solution to these equations with the minimal set of singularities 
in the physical strip 53m 6* G [0, vr]. It belongs to the class II in the classification of [1 



^m- r(i + a:)r(A-a;) r(|-x)r(i + A + x) 

^^^^-"r(i-x)r(A + a:)' '^^^^-r(i + x)r(i + A-x)- ^^^^^ 
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Figure 7: The two factorization equations for the backward scattering B. On each hne the 
scattering of an antiparticle (running downwards) with two particles (running upwards) de- 
composes into sequences of two-to-two scattering events controlled by the matrices S, F, B. 
Each line is associated to a rapidity that is carried either by a particle or an antiparticle 
according to the direction of the arrow. 
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In the above formulas x = iOjl-n and A = All other amplitudes may be obtained by 
using f ll29p and fll32p . This S-matrix has no singularity in the physical strip and thus the 
spinous do not form bound states. The overall sign, which is always left undetermined by 
the bootstrap program, was fixed by imposing that Siff) = — 1. This condition is expected 
to hold true for bosonic models in general [H]. Here we simply observe that this is the 
correct choice for matching with the large iV amplitudes of the x SU{N) models 

considered in this paper. 

The S-matrix fll35p has two remarkable properties. The first one concerns the spinon 
scattering phase in the symmetric channel S{6). It is easily seen to be identical up to an 
overall sign to the one found for the SU (N) chiral Gross- Neveu model if the substitution 
A —7- 1 — A is performed. We made use of this relation when computing the free energy 
density in the single-spin ground state, see Section. 14.21 The second remarkable feature is 
the scatteringphase of a fundamental excitation with spinon-antispinon pair sharing the 
same rapidityjj This phase is simply given by the product S{9)ti(9) and is identical to 
the scattering factor in the symmetric channel of the 0{N + 2) sigma model [12]. We did 
not make explicit use of this equivalence in this paper, but it would be relevant for the 
computation of the ground state energy density of the two equal spins solution mentioned 



at the end of the subsection 14.21 

Finally let us comment on the large rapidity expansion of S{6) and ti{6) 

S{e) r-. e-'^^ , ti(^)~e*^^, e>l. (136) 

This asymptotic behavior is indicating that the spinous behave as particles of fractional 
Lorentz spin, quite similarly to the case of the chiral Gross- Neveu model |5S]- To leading 
order at large N we have A — and the spinous are boson as expected. 



B.3 S-matrix with backward scattering 

We turn back to scattering with reflection when ri{6),r2{6) 7^ 0. To determine these extra 
factors we need to consider the Yang-Baxter equations for the two last processes in f ll32p . 
These are matrix equations that can be expanded over a basis of SU (N) covariant tensors. 
For generic values of N this leads to 5 + 5 independent functional relations for the S, F, B 
amplitudes. But not for N = 2. For this value of only 4 + 4 equations are available. 

^This state can be viewed as an adjoint bound state at threshold. 
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This is due to the following identity between Kronecker symbols, 



SikSjmSln — SinSjmSkl — ^^im^jn^kl — ^jk^lm^in + ^ik^lm^jn + ^^jk^ln^im , (137) 

which is valid only when the indices run through two values. This may also be understood 
from the viewpoint of group theory. For the processes considered the SU{N) decomposition 

.-r (N +2)N(N ~l) (N +l)N(N -2) , , 

N®N®N=^ — ® — 2 -®N®N, (138) 

applies only for N > 2. Each number on the r.h.s indicates the dimensions of the associated 
irreps and N/N in the l.h.s stand for the fundamental/antifundamental of SU{N). For 
N = 2 the irreps and N are isomorphic and the second channel in the r.h.s of fll38p is 
absent. It is natural to expect less constraints in this case. 

Let us now consider more closely these functional identities for the S, F, B amplitudes. 
Two of them are valid for all N 



ri(^^12)5(013)ri(^23) +tl(^12)ri(^13)tl(^23) = ^(^12)ri(^^l3)S(^23) , 
tl(^^12)5(^13)ri(^23) + ri(0i2)ri(^^i3)tl(^23) = 5(^^i2)tl(^^13)ri(^23) . 



(139) 



After renaming 

t,{e) = ST{e), r,{e) = SR{e), (i4o) 

they are revealed to be equivalent to the equations for the sine-Gordon theory This 
should not be surprising since the equations (I139p can be obtained by considering spinon- 
antispinon scattering in the U{1) subsector, which is spanned by Zi and Z]^. The solution 
to the U{1) factorized equations fll39p is known to be given by [12] 



.sinh(m-a6') .sinha6' 

S{9) = -I '-SR{e) , St{9) = I . ^ Sr{9) . 141 

sm a sm a 

The parameters a and a are arbitrary, though the latter is defined modulo 27r. This hints 
at the possibility of having solutions with continuous parameters. To verify this one has 
to also solve the equations accounting for the SU{N) degrees of freedom. This is where 
the aforementioned difference between N > 2 and N = 2 plays a crucial role. 

In the generic case N > 2, two of these additional equations are similar to fll39p with 
the only difference being the substitution S{9) = Ui{9) + Ui{9) — )• U2{9) — Ui{9). This 
however immediately excludes any continuous solutions since the ratio Ui{9)/u2{9) has 
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been already constrained in fll32p and is proportional to 9. The consequence is that the 
sine-Gordon type expressions fll4ip become rational and the continuous parameters are 
lost. In other words the extra SU{N) equations spoil the freedom in fll4ip . Assuming 
ri,r2,ti,Mi 7^ and taking into account the remaining SU{N) equations, one finds only 
two possible solutions 



corresponding to c = ±1. Unitarity and crossing yield then the additional relation 



The minimal solutions corresponding to c = ±1 are respectively given by the class III and 
IV solutions of [H]. The class III S-matrix is nothing else than the S-matrix for the 0{2N) 
sigma model, while the class IV, to our knowledge, has not yet emerged as a scattering 
matrix of an integrable sigma model. Thus, as pointed out by Berg et al. [H], there is no 
room for a one-parameter family of solution if > 2. In particular, for generic values of 

there is no interpolation between the refiectionless S-matrix and the 0(2A^) S-matrix, 
i.e. between type II and type III solutions in the classification of [14J. 

Interestingly enough, this conclusion is invalid for N = 2. Of course all solutions 
found at generic A^ exist also at A^ = 2. But with less number of constraints we expect 
a larger class of solutions, perhaps even with continuous parameters. This is in fact what 
happens and we have checked exphcitly that Wiegmann's S-matrix [30] solves all bootstrap 
equations for A^ = 2. The full solution reads 



cui{e) 



n{9) = U2{e) 



t2{0) = r2{9) , 



(142) 



S{9) 



Ssu{2){0)Sp{6) , 



(144) 



with the sine- Gordon soliton scattering phase 




(145) 



and the minimal SU{2) scattering phase 



r(i + f)r(| 
r(i-f)r(i 




(146) 
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All other amplitudes can be obtained by using fll4ip and fll32p . setting a = a = an, 
u = 71, and imposing crossing symmetry. For p = oo it is the 0(4) ~ SU{2) x SU{2) 
minimal solution while at p = 1 it becomes the minimal U{2) reflectionless S-matrix 
interpolating between class II and III of [H]. The overall sign in (11441) is chosen in accord 
with the generic p = 1 solution (11351) . 

We should mention here that there is a slight difference with the sine-Gordon S-matrix, 
in case of which the crossing symmetry requires that a = avr + vr in (I14ip . This is different 
from a = avr for = 2. The lack of the extra vr amounts to changing the sign of St{0) 
which in turn interchanges sectors which are even and odd under the action of the charge 
conjugation C. Mathematically, one has St^Sr — > StI^Sji. This modification is harmless 
as far as the physical properties of the model are concerned, but is important for attributing 
correct C-parity to the bound states in the N = 2 theory. 

C Counting argument : explicit calculation 

In Section 13. 2. II we have provided a synopsis of the counting argument proposed by Gold- 
schmidt and Witten [16]. We refer the reader to that paper for further details and expla- 
nations. Here we will apply this argument to the fermonic model. The following properties 
need to be observed when classifying the independent local operators. 

• Taking derivatives of the kinematical constraint zz = 1 and the Gauss law (12^ allow 
to discard the operators zD^z and their charge conjugates. 

• String of derivatives of opposite helicity can always be broken up. This follows from 
the equations of motion ( l38l) and the commutator 



• Similarly the helicity of the fermion should have the same sign as the derivatives 
acting on it. For example, the operators D^ip^ and their derivatives are forbidden. 

• The theory is invariant under charge conjugation. This property naturally extends 
to operators in both lists. 

Given the remarks above it is rather straightforward, if a little tedious, to compose lists A 
and B for the conservation law 



[D+, D_] ~ {D_zD+z - D+zD_z) . 



(147) 



d+{T__f = 0. 



(148) 
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We found that the set of anomahes consists of 16 elements 



A = ^DtzD+z + c.c. , {D'izD^z){D^zD+z) + c.c. , {D^zD^z){D+zD^z) + c.c. , 
{DlzD+z){D^zD^z) + c.c. , ip*_Diip_iplip+ + c.c. , D_ip*_D^_i;_ipli;+ + c.c. , 

, D^zD^zd-{i)*^i).)i)li)+ , 
i{D_zD+z - c.c.)D_zD^zip*_ip^ , i{DlzD+z + c.c.){ip*_D_tp_ - c.c.) , 
i{DlzD+z - c.c.)d_{ij*_ij^) , i{DlzD_z - c.c.)V';^+ , 
i{D_zD+z + c.c.){i)*_D'iip_ - c.c.) , i{D.zD+z - c.c.) {i^lDlip^ + c.c.) , 
i{D^zD+z - c.c.)D^i!*_D_ij_ , i{DtzD+z - c.c.)CV^-} , 

where c.c. denotes complex conjugation. The list of divergencies consists of two disjoint 
subsets 

B = B+UB^, (149) 

with 

B^ = d+ ^D^_zD_z + c.c. , D^zDlz , ii)*_Dti)_ + c.c. , iD_i)*_D'ii)_ + c.c. , 
iD_zD_z{^*_D_ij_ - c.c.) , i{D^_zD_z - c.c.)ip*_ij_ , 
{D_zD_zf , D_^*_D_i>_i)*_^i>^ I , 

B_ = d-^^DtzD+z + c.c. , D^zD^z{D^zD+z + c.c.) , {tp*_D^_^_ + c.c.)?/^;V+ 
D_ip*_D_i;_iplip+ , i{DlzD^z - c.c.)iplip+ , 
D_zD_z ij*_ilj_ipl'ilj+ , i{D+zD_z + c.c.){ij*_D_ij_ - c.c.) , 
i{D+zD_z - c.c.)d_{ip*_ip_) , i{DlzD+z - c.c.)ip*_ip_^ . 

It contains 17 elements. We should be careful with the interpretation of the results though. 
Since there are more operators in the set B than in the set A, a particular combination of 
the former operators should vanish at quantum level. This identity may be easily derived 
by differentiating the conservation law for the stress-energy tensor 

a+T__ + 9_T+_ = . (150) 

Note that the above equation, contrary to the one obtained by setting T+_ = 0, holds true 
also quantum-mechanically. Thus one element needs to be removed from the B-list. To 
check whether the conservation law fll48p may survive the quantization, we should also 
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strip off the combination 9+(T )^ or, equivalently, one of its constituent operators. This 

leaves us with 16 anomahes and only 15 divergencies, that is, with one unmatched anomaly. 
The above counting is generic and is oblivious to the fact that at special values of 
some operators may become redundant. This phenomenon might happen because our 
previous analysis does not guarantee that the kinematical constraint and Gauss law were 
exhausted. As it turns out, an exceptional relation among operators of the A-list may be 
found for N = 2, quite similarly to what we observed in Section lB.l.ll This relation can be 
understood as a fermionic deformation of the geometrical identity valid for the CP^ model 
[16]. Classically it is expressed as 

(D'^zD_z){D_zD^z) — [D'^_zD^z)[D_zD_z) + c.c. = fermionic operators, (151) 

where the r.h.s is a certain linear combination of (mixed) fermionic operators from the A- 
list, for example d^{D_zD^z)il)*_il)^il)*^il)j^, etc. Irrespectively of its precise form, which is 
likely to be quantum corrected anyway, such a relation allows us to eliminate one operator 
from the A-list. This is the sole modification and the B-list is valid as it stands. One 
thus infers that a non-trivial higher-spin conserved current should exist for N = 2. This 
mechanism was first observed for the CP^ model [IB]. Can we find relation similar to fll5ip 
for N > 21 We believe this is not possible for the following reason. If such a relation existed 
it should be a fermionic deformation of an identity from the CP^~^ model. The latter 
however does not exist [IB] and we expect the above counting to be correct for any N > 2. 

Finally we stress that the number of unmatched anomalies is independent of the value 
of the Thirring coupling and the fermionic charge. This number also coincides with its 
counterpart for the CP^~^ sigma model The important difference however is that for 
the fermionic extension we have one continuous parameter at our disposal. Can one pick 
its value such that the coefficient of the unmatched anomaly is cancelled for N > 21 The 
answer seems to be positive as the results in other sections of this paper indicate. 
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